Three fully polarized fermions close to a p-wave Feshbach resonance 
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We study the three-body problem for three atomic fermions, in the same spin state, experiencing 
a resonant interaction in the p-wave channel via a Feshbach resonance represented by a two-channel 
model. The rate of inelastic processes due to recombination to deeply bound dimers is then esti- 
mated from the three-body solution using a simple prescription. We obtain numerical and analytical 
predictions for most of the experimentally relevant quantities that can be extracted from the three- 
body solution: the existence of weakly bound trimers and their lifetime, the low-energy elastic and 
inelastic scattering properties of an atom on a weakly bound dimer (including the atom-dimer scat- 
tering length and scattering volume) , and the recombination rates for three colliding atoms towards 
weakly bound and deeply bound dimers. The effect of "background" non-resonant interactions in the 
open channel of the two-channel model is also calculated and allows to determine which three-body 
quantities are 'universal' and which on the contrary depend on the details of the model. 
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I. INTRODUCTION 

Fcrmionic superfluidity with p-wavc pairing is related 
to a large class of subjects in very different areas of 
physics including condensed matter, astrophysics and 
particle physics [1, 2]. As already observed in '^He exper- 
iments, the phase diagram in these systems can be very 
rich [3]. Moreover, the possible observation of quantum 
phase transitions together with the existence of exotic 
topological defects in p-wave superfluids bring a lot of 
interest in their study. 

Presently, there is some hope that p-wave superfluidity 
and its intriguing properties can be observed with ultra 
cold atoms [4]. Indeed, thanks to the concept of Feshbach 
resonance [5] , it is possible to tune the inter-atomic inter- 
action and to achieve strongly correlated regimes in ul- 
tra cold dilute atomic gases. First realized in the s-wave 
channel with bosonic species [6, 7], the Feshbach reso- 
nance is currently used for achieving BEC-BCS crossover 
experiments for the two-component Fermi gas in a regime 
of temperatures where the system can be superfluid [8- 
23] . In one spin component Fermi gases, as a consequence 
of the Pauli exclusion principle, two-body scattering pro- 
cesses are forbidden in the s-wave channel and at low 
temperatures are dominant in the p-wave channel. The 
two-body cross-section which is usually negligible in this 
channel can be greatly enhanced using a p-wave Feshbach 
resonance. This resonant regime is now obtained for ^''K 
[24-27] and ^Li atoms [28-30]. The production of p-wave 
shallow dimers in ultra cold ^Li [28] and gases [27] 
opens very interesting perspectives for the realization of 



a superfluid p-wave phase. 

In these experiments, an external magnetic field tunes 
the energy of a two-body p-wave bound state in a closed 
channel and for a small detuning with respect to the open 
channel, a resonance occurs in a two-body p-wave scat- 
tering process. Moreover, due to the presence of a mag- 
netic field, the interaction strength depends on the orbital 
channels considered -a major difference with respect to 
what happens in superfluid "^He [25]. As a consequence, 
the question of the symmetry of the low temperature 
ground state in one component fermionic species is non- 
trivial. Studies of this many-body problem are essentially 
mean-field and depending on the experimental realiza- 
tions, they predict the occurrence of p^: -I- ipy (axial), Px 
(polar) or intermediate phases [31-36]. These predictions 
lead to possible studies of quantum phase transitions in 
such systems. However, the main issue in the achieve- 
ment of a p-wave superfluid concerns atom losses which 
arc large in present experiments [24-30]. The question 
whether or not it is possible that the system thermalizes 
is then a crucial point. 

Concerning s-wave resonant Fermi systems, few-body 
studies have proven to be very successful in understand- 
ing properties of the superfluid gas in the BEC-BCS 
cross-over [37-41]. These studies explain the large life- 
time of the system observed at resonance and also predict 
the dimer-dimer scattering length which is involved in the 
equation of state for the dilute BEC phase. Surprisingly, 
although general many-body properties are rather well 
known in p-wave superfluids -thanks to contributions 
from the condensed matter community, few body prop- 
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erties in these systems have been less studied [42, 43]. 
However, following the example of the works done in the 
s-wave channel, few-body problems for p-wave pairwisc 
potential are valuable for a determination of properties 
in the strongly interacting dilute gas beyond a mean-field 
analysis. As an example, we note that consequences of 
the existence of trimers first found in the present work 
have already been taken into account for an estimation 
of the lifetime of p-wavc shallow dimers [44] . 

In this paper, we consider three identical fermions close 
to a p-wave Feshbach resonance. We determine their low 
energy scattering properties together with the possible 
existence of trimers. Our study is also a first step to- 
ward an understanding of the atom losses observed in 
present experiments [24, 26-29]. The paper is organized 
as follows. In section II, we recall basic properties, for an 
isotropic short range interaction, of resonant two-body 
p-wave scattering processes [45] . In the resonant regime, 
two parameters arc needed for a description of the low en- 
ergy two-body properties: the scattering volume and 
also the p-wave equivalent of the effective range parame- 
ter hereafter denoted by a. For large and positive values 
of Vs there exists a shallow p-wave dimer of internal an- 
gular momentum one, that is with three-fold degeneracy. 
For a potential with a compact support of radius 6, we 
show that at resonance (Vs = oo) the effective range pa- 
rameter cannot reach arbitrarily small values and ab > I. 
Consequently, unlike what happens in s-wavc resonances, 
there is no scale invariancc at low energy and a unitary 
regime cannot be obtained via a p-wave resonance [46] . 
In section III, we introduce the main model Hamiltonian 
that we use in this work. It is a two channel model of the 
p-wave Feshbach resonance [30] where free atoms in the 
open channel interact with a molecular p-wave state in 
the closed channel, of threefold degeneracy provided that 
one neglects the effect of the dipole-dipole interaction 
in presence of the Feshbach magnetic field. The inter- 
channel coupling amplitude, as a function of the relative 
distance of the two atoms, is a Gaussian of range b which 
mimics the van der Waals range of a more realistic two- 
body potential. We first briefly determine the two-body 
coUisional properties of this model. At large coupling the 
resonance is broad afo ~ 1 and for Vs large and positive 
the shallow dimer is essentially in the open channel. In 
the opposite regime for a weak coupling, the resonance is 
narrow, ab ^ 1 and the shallow dimer is almost entirely 
in the closed channel. In section IV, we derive an inte- 
gral equation for the three body problem. We consider 
solutions of total angular momentum J = 1 and by using 
the rotational symmetry of the Hamiltonian, wc reduce 
the problem in each involved symmetry sector (odd or 
even) to a one dimensional integral equation. In both 
sectors, we predict the existence of one trimer for suffi- 
ciently broad resonances. These trimers can exist in a 
regime where there is no shallow dimer (for large and 
negative values of the scattering volume) and are inter- 
esting examples of Borromean states [47] , since we find 
that they arc not linked to an Efiniov effect, contrarily to 



[43] . We determine also the atom-dimer scattering length 
aad as a function of the effective range parameter a and 
the potential range b, for different values of the scattering 
volume. At resonance (Vs = oo), aad takes large values 
(that is significantly larger than the potential range b) 
only in the vicinity of the threshold of existence of a 
trimer. The recombination rate of three incoming atoms 
into a shallow dimer and one outgoing atom is computed; 
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it is shown analytically to vary as Vs for large values 
of the scattering volume, away from the trimer forma- 
tion threshold; this differs from the Vs^'^ law put forward 
in [42] on the basis of a dimensional analysis ignoring a 
possible contribution of the effective range parameter a, 
but is still compatible with the numerics of [42]; finally, 
the recombination rate is shown analytically to present a 
Fano profile as a function of a close to this trimer thresh- 
old. In section V, we calculate the losses due to the re- 
combination into deeply bound dimers. Since these losses 
are not present in our model Hamiltonian, we estimate 
them from the probability that three atoms are within 
a volume of the order of b'^ and we obtain the lifetime 
of trimer states, the losses due to atom-dimer inelastic 
scattering and the three body recombination rate toward 
deep molecular states from asymptotically free atoms. 
Finally we make the model more realistic by including 
an attractive interaction in the open channel, in addi- 
tion to the coupling with the closed channel, in section 
VI: we recalculate the trimer energies, the atom-dimer 
scattering length and the recombination rate to weakly 
bound dimers, and we physically explain the impact on 
these quantities of a non-resonant interaction in the open 
channel. Wc conclude in section VII. 

II. BASIC PROPERTIES OF THE TWO-BODY 
P-WAVE SCATTERING 

A. The scattering amplitude 

We consider in this section two particles of mass m 
in the same spin state and in the center of mass frame, 
scattering in free space via a rotationally invariant short 
range interaction potential. We assume for simplicity 
that the interaction potential scatters only in the p-wave 
channel, so that at large distances, where the effect of 
the potential is negligible, the scattering wavefunction of 
energy E = h^k'^/m, fc > 0, takes the form 

^.(r).e- + 3/(fc)k.4(^) (1) 

where r is the relative position of the two particles, ±k 
are their incoming wave-vectors, and we have introduced 
the unit vectors f = r/r and k = k/fc. The function f{k) 
is the so-called reduced scattering amplitude since the 
angular dependence of the scattered wave has been pulled 
out. We note that Eq.(l) becomes exact (that is one 
can replace ~ by =) for a compact support interaction 
potential, when r is out of the support of the potential. 
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In this subsection, we briefly review some basic proper- 
ties of this p-wave scattering ampUtude /(fc). As a con- 
sequence of the unitarity of the S'-matrix of scattering 
theory, it obeys the optical theorem, 

lmf{k)^k\f(k)\^ 



which imphcs 



fik) 



1 



(2) 



(3) 



u{k) + ik 

where u{k) is a real function. For cold atoms, the low 
energy scattering properties are crucial and we assume 
that u{k) has the following low- A; series expansion, 

uik) = ^^+a + 0{k^). (4) 

The so-called scattering volume plays a role similar to 
the scattering length in the s-wave channel: the resonant 
situation corresponds to the limit \Vs\ — > oo. 

Another crucial property of the reduced scattering am- 
plitude is that its analytic continuation to negative ener- 
gies, that is to imaginary values of k, gives information 
on possible bound states in the two-body problem, in 
the form of poles of f{k). More precisely, setting k = iq, 
where q > 0, the solutions (?dim > of the equation 

correspond to bound states of the scattering potential, 
that is here to dimers of rotational quantum number S = 
1, with a binding energy 

Edim = h^qli^/m. (6) 

The wavefunction of such a dimer, "out" of the potential 
(again, this has an exact meaning for a compact support 
potential), is a solution of the free Schrodinger's equation 
in the p-wave channel, so that we may take it of the form 

1/2 



47r j r dr 



(7) 



where A/" is a normalization factor and is the compo- 
nent of r along direction 7 = cc, y or z [48]. 

The knowledge of the dimer wavefunction "inside" the 
potential requires a full solution of Schrodinger's equa- 
tion. However it is possible to access the normalization 
factor M directly from the knowledge of the scattering 
amplitude. Using the closure relation 

/ (03V'k(rM(r') = 5(r-r')-^0.(r)0*(r'), (8) 



in the limit of large r and r' 

w? = - 



we obtain 
2 



^ (9) 

9dim [1 - jw'(jgdim)] ' 

assuming that u{k) has a series expansion with even pow- 
ers of k only and using contour integration in the complex 
plane to single out the contribution of the poles of f{k). 
As we shall see, this relation (9) may be used to put 
constraints on the parameter a. 



B. Constraint on the parameter a close to 
resonance 

Whereas the scattering volume can be adjusted at will 
by a Feshbach resonance driven by a magnetic field, the 
value of a on resonance cannot be adjusted the same 
way so it is important to determine what are its possible 
values on resonance. 

We assume that ofj-cs 7^ 0, where arcs is the value of a 
on resonance [49]. In the resonant limit, we see from the 
low-fc expansion oiuik) that there exists a weakly bound 
dimer on the side a^cs^s > of the resonance [50]: 



1 



9dii 



V arcs Vs 



From Eq.(4) and Eq.(9) we obtain 



(10) 



(11) 



This imposes ares > 0. This is in sharp contrast with 
the case of s-wave scattering, where the effective range 
Te can take any sign on resonance. 

For a compact support potential, vanishing outside a 
sphere of radius 6, that is for r > 6, the normalization of 
the dimer wavefunction to unity imposes 



dV|0(r)|2 < 1. 



(12) 



r>b 



Calculating the resulting integral with the expression 
Eq.(7) leads to [51] 



1 



< 1. 



In the limit \Vs\ -^00 this leads to [52] 



1 

> - 
- b 



(13) 



(14) 



where, again, arcs is the value of a on resonance |Vs| = 
00. In the zero range limit 6—^0, we see that arcs cannot 
tend to zero, but on the contrary has to diverge! This 
is in sharp contrast with the s-wave case, where one can 
find models for the interaction potential where — > in 
the zero range limit 6^0. 

To illustrate these properties on a simple example, we 
give in Fig. 1 the values of a and for a square well inter- 
action potential, as functions of the well depth. We see 
on the figure that (14) is satisfied at resonance, and that 
a is no longer constrained by this condition away from 
resonance, and may even vanish and become negative. 



III. MODELING OF THE RESONANT p-WAVE 
INTERACTION 



In this section, we introduce the main model used in 
this paper to describe the p-wave interaction between 
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FIG. 1: For a square well interaction potential V{r) — 

^ 9{b — r), where 9 is the Heaviside function, values of a 

in units of 1/6 (dashed line) and Vs in units of (solid line) 
as functions of kob. Note the divergence of a when Vs = 0. 



same spin state fermions close to a resonance. It is sim- 
ply a two-channel model of a Feshbach resonance, that is 
a direct generalization of the s-wave two-channel model 
[53] to the p-wave case, in the spirit of [30] . It is extended 
in section VI to include direct interactions among atoms 
in the open channel. 



Hamiltonian written in second quantized form: 



H = 



A 



(2^ 

d^kcPk 
{2nf 



21.2 



-7; Ok^k 

2m ^ 



mol 



Am 



k-k' 



h.C. 



(15) 



The annihilation and creation operators for fermions 
(that is for the atoms in the open channel) in plane waves 
of wave vectors k and k' obey the anticommutation rela- 
tion 



{ak,4'} = (2^)''5(k-k') 



(16) 



which corresponds to the convention (r|k) = e*'^ '" for the 
plane-wave. The operator 6.^^k annihilates a molecule (in 
the closed channel), with a center of mass momentum 
?ik, in one of the three degenerate internal states 7 in 
the Snioi = 1 molecular rotational manifold; we take here 
for 7 one of the directions x, y or z, which amounts to 
using the chemistry basis {[7)}, where [7) is an cigen- 
statc of zero angular momentum along direction 7, rather 
than the standard basis {|m = 0, ±1)}- As we mentioned, 
molecules are treated as bosons so that the 6's obey com- 
mutation relations 



Vk,6y,k'] - V(27r)'<5(k-k'). 



(17) 



A. Model Hamiltonian 



As is standard in a two-channel model, the atoms may 
populate either the open channel, where they are treated 
explicitly as fermionic particles, or the closed channel, 
where they exist only under the form of specific tight two- 
body bound states, here referred to as molecules; these 
molecules are treated as bosons, and have an internal ro- 
tational state of spin 511101 = 1 since they are p-wave two- 
body bound states. We assume that the three rotational 
sublevcls of a molecule are degenerate: even if this is not 
exactly true in practice because of the effect of the dipole- 
dipole interaction in presence of the magnetic field used 
to produce the Feshbach resonance [25], this will make 
our model rotationally invariant and greatly simplify the 
algebra for the three-body problem. For simplicity, we 
also assume that there is no direct interaction among 
the fermionic particles, the resonant p-wave atomic in- 
teraction being taken into account through the coupling 
between fermions and molecules. As already mentioned, 
this simplifying assumption is removed in section VI. 

The situation is represented schematically in Fig. 2. 
Mathematically, it corresponds to the following free space 



Also, the a, fermionic operators commute with the 
bosonic ones b and . In addition to its center of mass ki- 
netic energy, each molecule has an internal energy E'moi, 
defined in the absence of coupling between the open and 
the closed channels, and counted with respect to the dis- 
sociation energy of the open channel. 

Whereas the first contribution in the right-hand side 
of (15) simply corresponds to non-interacting gases of 
atoms and molecules, the second contribution describes 
the coupling between the two species, that is between 
the open and closed channels, responsible for the p-wave 
resonance. This inter-channel coupling depends on the 
relative momentum between two atoms through the func- 
tions x-y'i here, we are in the case of a p-wave coupling so 
we take 

X(k)=ke-'=^^V2 (18) 

where b is the range in real space of the inter-channel 
coupling, of the order of the radius of the closed chan- 
nel molecule. The overall amplitude of the inter-channel 
coupling is measured by the coupling constant A, taken 
here to be real; it has not the dimension of an energy, 
but rather has the same dimension as fi^b^/'^/m. As we 
already mentioned, the model is summarized in Fig. 2. It 
holds at low kinetic energies, below the dissociation limit 
Voo of the closed channel. 
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FIG. 2: Schematic view of a Feshbach resonance configura- 
tion: the atoms interact via two potential curves, plotted as a 
function of the interatomic distance. Solid line: open channel 
potential curve. Dashed line: closed channel potential curve. 
When one neglects the coupling A between the two curves, the 
closed channel has a molecular state of energy -Emoi with re- 
spect to the dissociation limit of the open channel. Note that 
the energy dependence of the two curves is purely indicative, 
and the spacing between the solid curve and the dashed curve 
was greatly exaggerated for clarity. 



B. Two-body aspects 

Before solvmg the three-body problem, it is important 
to understand the two-body aspects of the model, in the 
form of the reduced scattering amplitude f{k) and the 
related properties of possible dimcrs, according to the 
general discussion of section II. 

We thus calculate the scattering state of two atoms in 
the center of mass frame, that is for a zero total momen- 
tum. The most general state vector is thus a coherent 
superposition of two atoms (in the open channel) and 
one molecule (in the closed channel): 



(2^ 



A{k)aial^\0). (19) 



Since the molecule has a zero total momentum, its state is 
characterized by the three complex amplitudes P-y in each 
of the internal rotational states "/ = x,y and z. On the 
contrary, the two atoms can have opposite but arbitrary 
momenta k and — k, hence the a priori unknown function 

Injecting this ansatz in Schrodingcr's equation (E — 
H)\'i>) = 0, and projecting onto the molecular subspace 
and the atomic subspace respectively, one finds that 
Schrodingcr's equation is satisfied when A and (3 satisfy 



(20) 



E 



A(k) + A/3 • x(k) = (21) 



Equation (21) does not specify A(k) in a unique way, 



for a positive energy E, since E — h'^k'^/m vanishes for 
some values of k. To obtain the scattering state of two 
atoms, one takes a more specific form of the ansatz, cor- 
responding to the superposition in the open channel of an 
incoming wave of wavevector ko and a purely outgoing 
scattered wave, 



A(k) = (27r)3,5(k - ko) - A 



/3 • X(k) 



E + iO+ - 



(22) 



Here i5 = ?i^fcg/m > is the total energy of the scatter- 
ing state. 

The general scattering theory [54] relates the scattering 
state to the incoming state |*o) by |*) = (1 + Gor)|*o) 
where T is the T-matrix and Gq the resolvent of the 
non-interacting Hamiltonian. From this identity it is 
then apparent that the matrix element of the T-matrix 
in Fourier space is related to the numerator of the last 
term of Eq.(22): 



(k|r(i? + zO+)|ko) = -A/3.x(k). 



(23) 



From the known relation between the scattering ampli- 
tude and the T-matrix [54] , we get the reduced scattering 
amplitude 



I{k^) 



3{E-E„ 
2A^ 



E+iQ+- 



k^e-''^^ 



(24) 



The choice of the Gaussian envelope in x(k) allows an 
explicit expression for the scattering amplitude. After 
complexification of fco by analytic continuation, setting 
^0 = ^^O; 9o > 0, we obtain 



1 



An 



+ 



2m2A2 
d^k k'^e-''^^ 



{ql + mT'moi/?!^) 



(2^)3 ql- 



1 + 



fc2 



(25) 

9ocrfc (go^) (26) 



where erfc is the complementary error function that tends 
to unity in zero. With this complexification technique, 
it is straightforward to identify the parameters Vs and 
a appearing in the low-A: expansion (4) and to get the 
explicit expressions: 



1 


1 


Vs ^ 


27ri/263 






a = 





E„ 



7ri/26 



2^2- 



(27) 
(28) 
(29) 



This illustrates the fact that one can tune Vg to — oo or 
+00 by shifting the molecular energy -Emoi (in practice 
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with a magnetic field B) around the value E^^^ such 
that the right hand side of (27) vanishes, i?moi — ^moi — 
KB -Bo). 

We have introduced the convenient quantity arcs, 
which is the value of a exactly on the Feshbach reso- 
nance. We see that arcs depends on the inter-channel 
coupling A, and is bounded from below by the inverse of 
the potential range, within a numerical factor depending 
on the details of the model, here I/tt^/^ . In principle, arcs 
can take any possible value above this limit, depending 
on the value of the interchannel coupling A; in practice, 
of course, A is not easily tunable so arcs is fixed for a 
given experimental configuration. 

By a direct generalization of a well established s-wave 
terminology, we may classify the p-wave Feshbach reso- 
nances as 

• a broad resonance {A:$>h'^b^/^/m): arcs ~ 1/6 

• a narrow resonance (A < h^b^^^ /m): arcs ^ 1/6. 

We recall that this terminology can be motivated as 
follows: If one assumes that i^rnoi is an affinc func- 
tion of the magnetic field B with a slope /i, and that 
= Vs^[l - AB/{B - Bo)] in a more complete theory 
including the fact that takes a finite value V^^ far from 
the resonance (due to the direct interaction in the open 
channel, neglected here) and presumably of the order of 
6'^, one finds a resonance width 



(30) 



It remains to compare this resonance width to the 'natu- 
ral' energy scale h^/mb"^ to obtain the above mentioned 
terminology. 

The last point to discuss for the two-body problem is 
the existence or not of a two-body bound state in the 
open channel. We shall refer to such a bound state as a 
dimcr, in order not to confuse it with the molecular state 
in the closed channel. Mathematically, such a dimer is a 
zero of 1/ f{iqo) with go > 0. The expression in between 
square brackets in the right hand side Eq.(25) is a de- 
creasing function of qq that tends to ~oo for go ~* -l-oo. 
Hence there exists at most one dimer in our model Hamil- 
tonian. There exists one if and only if the expression be- 
tween square brackets is positive in go = 0, that is if and 
only if >0 [55]. 

When a dimer is present, one can express analytically 
its wavcfunction 0(r) in the open channel, in terms of 
exponential and crfc functions, and one can calculate the 
occupation probability of the closed channel, Pcioscd — 
|/3p after proper normalization of \'i>) in the center of 
mass frame [56]; 



(pk 

(2^ 



|A(k)|2 = 1. 



(31) 



An equivalent way to obtain |/3|^ is to calculate the large 
r behavior of 0(r), which is proportional to 13 and which 



is related to the normalization factor M in Eq.(7), and 
then to use the general relation Eq.(9) [57]. Both ways 
lead to the expression 



1 TO^A^ e'Jdi. 



Pcloscd Ml^ 



(32) 



The value of Pciosed for an infinite scattering volume can 
be cast in the very simple forms 



7) 

inclosed 



m^A^ 
1 - 



1 



(33) 
(34) 



The expression (33) is quite remarkable since it is 'uni- 
versal': It does not involve the interaction range b and, 
as we have checked, it is not specific to the choice of a 
Gaussian cut-off function in x(k). It was already derived 
in [30], see the unnumbered equation following equation 
(9) of that reference. In the vicinity of the resonance, we 
see on the expression (34) that, in the dimer wavcfunc- 
tion, the closed channel is strongly occupied for a narrow 
resonance and is weakly occupied for a broad resonance; 
Pcioscd tends to zero in the broad resonance limit. 

To conclude this review of the two-body aspects, we 
point out a striking property of the dimer, very different 
from the usual s-wave case: In the limit Vs/b^ +oo, 
we find that the dimer wavefunction (j>{r) has a well de- 
fined, non-zero limit, tending to zero as 0(l/r^) at large 
r. This can be directly seen in momentum space: for 
gdim = 0+, the function A(k) is 0{l/k) at low fc, which 
is indeed square integrable around the origin k = 0. In 
other words, at the threshold for the formation of the 
dimcr, the dimer wavefunction is a well defined non-zero 
and square integrable function. 



IV. SOLUTION OF THE THREE-BODY 
PROBLEM 



This is the central section, where we solve the three- 
body problem within the two-channel model close to a 
p-wave resonance. The mathematical structure of the 
model, with a single molecular state occupied in the 
closed channel and no interaction potential in the open 
channel, is such that the three-body problem is amenable 
to an integral equation for a one-body 'wavefunction'. 
This integral equation becomes easily solvable numeri- 
cally if one further uses the rotational symmetry of the 
Hamiltonian. We then obtain predictions for three phys- 
ical situations, (i) the existence of three-body bound 
states, that is of trimers, (ii) the scattering of an atom on 
a dimer, and (iii) the scattering of three atoms, leading 
to recombination processes, that is to the formation of a 
weakly bound dimcr and a free atom. 
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A. Derivation of an integral equation 



This equation is satisfied for the choice 



We start with the most general ansatz for the three- 
body problem in the center of mass frame, that is for a 
zero total momentum. Because of the conversion of pairs 
of atoms into molecules and vice-versa, the three-body 
ansatz is a coherent superposition of three fermions (all 
three atoms in the open channel) and of one molecule 
plus one fermion (one atom in the open channel and two 
atoms tightly bound in a molecule in the closed channel) : 



I*) = 



l3-,{K)bi^al^\Q) 



AiK,k)a\^^^a\^_^al^\0).i35) 



The one molecule plus one fermion part is parameterized 
by three one-body 'wavefunctions' /3^, here in Fourier 
space; wc shall derive an integral equation for them. The 
three fermion part A can be parameterized by two Jacobi- 
like coordinates in momentum space since the total mo- 
mentum is zero. For pure convenience, we impose that 
A(K, k) is an odd function of k, to reduce the number of 
terms involving A in the integral equation for (3. 

We inject the general ansatz for l^*) in Schrodingcr's 
equation {E — H)\'i') — 0, where the total energy E is 
at this stage of arbitrary sign. Projecting Schrodinger's 
equation on the subspace with one molecule and one 
fermion gives an equation for f3 with the function A ap- 
pearing in a source term: 



E- E, 



niol 



4m 



A{K,k) + 2A 



rK + k, 



2A 
3 



K 



d^k 
(2^ 
1 



X*(k) 
= 0. (36) 



Projecting Schrodinger's equation on the subspace with 
three fermions leads to: 



d^Kd^k 
(27r)6 

+A/3(K).x(k)|al„..a',,, .a 



K+k"iK-k"-K 



A(K,k) 

= 0. (37) 



A(K,k) = Ao(K,k) 



A/3(K) • x(k) 



E + iQ+ 



(38) 



For a positive total energy > 0: ^0 represents a possi- 
ble incoming wave of three free atoms, and it is an eigen- 
state of the kinetic energy operator in the center of mass 
frame with energy E\ in presence of such an incoming 
free wave, the second term in A represents the scattered 
wave in the open channel, which is guaranteed to be out- 
going by the standard substitution E ^ E + iO+. As 
we have imposed the convention that A{'K, k) should be 
an odd function of k, one has to apply the same conven- 
tion to ^o(K,k); note that the last term of Eq.(38) is 
automatically an odd function of k, since x(k) is. For 
a negative total energy E < 0, the expression between 
square brackets in (37) can not vanish, Aq = and the 
-f-iO"'' in the denominator can be omitted. 

Injecting Eq.(38) in Eq.(36), wc obtain an integral 
equation for (3. The term j4(K,k) of Eq.(36) gives a 
contribution simply proportional to /3(K), with a in- 
dependent factor; collecting it with the factor in between 
square brackets in the first term of Eq.(36) gives a K 
dependent expression that can be recognized as being 
proportional (with a K dependent factor) to the inverse 
of the scattering amplitude of two atoms at the energy 



E, 



rcl 



E 



4m 



(39) 



with the determination fcici > for E^^i > and k-^d/i > 
for £"101 < 0. This relation can be seen as a consequence 
of the Jacobi-like parameterization of the momenta of the 
three fermions that we have used in (35): if three free 
fermions of total energy E have momenta ikid + K/2 
and — K, then the modulus fcroi will obey (39). 

We finally obtain the general integral equation for the 
/3^(K): 



3/(fcrel) 



/3(K) + 87r 



d^k , 



\. I A /3(k)-x(K + lk) 

2 J + k^ +K-k-'m{E + iO+)/h^ 



mA 



d^k 
(2^ 



X*(k) 



^o(K, k) + 2Ao ( -^K + k, -^K - Ik 



(40) 



r 



In what follows we shall solve this integral equation for energy (if there exits a dimer); then Aq = and /3(K) 
various physical situations, (i) In the search for trimcrs, is not subjected to any specific boundary condition, (ii) 
one assumes an energy E below zero and below the dimer In the low energy scattering of an atom on a dimer, the 
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energy is above the dimer energy but still negative; then 
Aq = and one has to introduce a specific ansatz for 
/3(K) to enforce the boundary conditions corresponding 
to such a scattering experiment, (iii) In the scattering 
of three incoming atoms, the total energy is now non- 
negative so that ^0 7^ 0; we shall assume that this scat- 
tering experiment is performed for > so that there 
exists a dimer in the two-body problem, which can be 
formed by a recombination event in the three-body scat- 
tering; then one introduces an ansatz for /3(K) describing 
the presence of a purely outgoing wave of such a dimer 
(with an opposite momentum atom). 



which corresponds to the even state with J = 1 and 
vanishing angular momentum component along x axis. 
After some calculations Eq.(40) can be turned into an 
integral equation for i?L=i (in the even sector) or into 
coupled integral equations for i?L=o and i?L=2 (in the odd 
sector), as detailed in the Appendix A. The remaining 
unknown functions depend on a single real variable K so 
that a numerical solution is reasonable. 



C. Existence of weakly bound trimers 



B. Symmetry sectors from rotational and parity 
invariance 

Formally Eq.(40) is an equation for a spinor /3(K), 
with an internal spin 5*1001 = 1; here this internal spin 
corresponds to the rotational degrees of freedom of the 
molecule (in the closed channel); the orbital variable K 
here corresponds to the relative atom-moleculc momen- 
tum. The homogeneous part of Eq.(40) is invariant by a 
simultaneous rotation of the spin and orbital variables of 
the spinor. The total momentum J, obtained by addi- 
tion of the spin Smoi and the orbital angular momentum 
L, is therefore a good quantum number. In this paper, 
we shall restrict to the manifold J = 1 , which can be ob- 
tained from L ^ plus S'moi = 1, or L = 1 plus 5moi = 1, 
or L = 2 plus S'moi = 1- In addition, the homogeneous 
part of Eq.(40) is invariant by parity (combining the par- 
ity on the spin variables and on the orbital variables). 
This decouples the J = 1 manifold in two sectors, 

• even sector: L = 1 plus 5moi = 1 

• odd sector: L ~ plus 5moi = 1 and L = 2 plus 

S'moi — 1 

Applying the standard algebra of addition of angular 
momenta, we obtain the following ansatz in the odd sec- 
tor. 



/3(K) = BL=o{K)e. 



Bl=,{K)^^K, (41) 



where is the unit vector along z axis. This ansatz cor- 
responds to a total angular momentum J ~ 1 with van- 
ishing angular momentum component along z, mj = 0. 
Considering the other components mj = ±1, or equiva- 
lently the states with vanishing angular momentum com- 
ponent along X and along y respectively, would lead to 
equivalent results, as guaranteed by the rotational invari- 
ance of the Hamiltonian. 

Similarly, we take as ansatz in the even sector 

;3(K) = :?^^KAe. 



K 

Bl=i{K) 



[(K . e,)e, - (K • e,)e,] (42) 



We investigate here the existence of three-body bound 
states, that is of trimers, in our model Hamiltonian. 
These trimers have of course a negative total energy E. If 
one is on the Vs > side of the resonance, where a dimer 
of energy — i^dim exists, one further has E < — i?dim to 
have stability of the trimers with respect to dissociation 
into an atom and a dimer; if this condition was not satis- 
fied, the trimers would not exist as true stationary states 
but would rather be resonances in the atom-dimer scat- 
tering process. 

These constraints on the energy have the following 
mathematical consequences. Since E < 0, the source 
term Aq in (38) is identically zero, so that Eq.(40) be- 
comes homogeneous. Since E < — Sdim, on the side 
Vs > of the resonance, the scattering amplitude f{k^ei) 
in the denominator of the first term of Eq.(40) is non 
zero for all K and the linear operator L{E) represent- 
ing the integral equation has a smooth action over the 
class of regular /3(K) functions. Numerically, one can 
then discretize the variable K with no particular care, 
and approximate L{E) by a matrix. The existence of 
a trimer corresponds to a non-zero-dimension kernel of 
the operator L{E); in practice, we look for the values of 
E such that the approximating matrix has a vanishing 
eigenvalue. The explicit form of L{E) for the ansatz in 
the even and odd sector can be deduced from the ap- 
pendix A. In the same appendix, it is also explained how 
to normalize the state vector of the trimer. 

For values of |Vs| ^ 6"^, we have found either zero or 
one trimer in each symmetry sector (with threefold ro- 
tational degeneracy when the trimer exists) . The energy 
of the trimer is written as —h?q^^i^/m. Then gtrim as a 
function of ab is given in Fig. 3, for the even and the odd 
sectors. We found no evidence of Efimov effect: in each 
symmetry sector, we found at most one trimer, and there 
is no oscillation of the /3(K) with K a,s a. function of K, 
see Fig. 4. 

We note that, in real experiments with atoms, these 
trimers may acquire a finite lifetime, due to the formation 
of deeply bound dimers by three-body collisions. This 
process is not contained in our Hamiltonian, since H does 
not support deeply bound dimers for |Vs| 3> 6^; its rate 
is estimated by a simple recipe in subsection VA. 
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0.5 1 




0.75 



a b 



FIG. 3: (Color online) For fixed values of the scattering vol- 
ume Va, parameter gtrim of the trimer (when it exists) as a 
function of ab; gtrim is related to the negative energy — -Etrim 
of the trimer by i^trim ~ ^^'j'?rim/"^. (a) Even sector, (b) 
odd sector, as defined in subsection IV B. Solid line (black): 
|Vs|/6^ — oo. Above the solid line, positive values of Vs- 
short dashed line (blue): Vs = 10* fe"'; dashed line (red): 
Vs = lO^fe^; dashed-dotted (green): = 1006^. Below the 
solid line, negative values of Vs'- short dashed line (light blue): 
Vs = -W^b^; dashed line (orange): = -lO^b^; dashed- 
dotted (dark green): Vs = — 100&^. At the threshold for the 
existence of the trimer as a true bound state, on the Vs > 
side of the resonance, where a dimer exists, the trimer bind- 
ing energy vanishes, so that the energy of the trimer coincides 
with the one of the dimer, and qtrim = Qdim (see text). 

D. Atom-dimer scattering 

We consider here the scattering problem of an atom on 
a dimer, which corresponds to the positive Vs side of the 
resonance and to a total energy E > — i^dim- For simphc- 
ity, we restrict to the low energy limit of this scattering, 
with a relative kinetic energy of the incoming atom and 
the dimer much smaller than the binding energy of the 
dimer: 

i^ + £^dim<-Bdim. (43) 

As a consequence, the total energy is negative, so that 
energy conservation prevents the dimer from being dis- 



(a) 




3 4 



K[l/b] 




(b) _ 



1 2 3 4 5 

K [1/b] 

FIG. 4: For Vs = oo, and ares = Qth [right on the thresholds 
for the existence of a trimer, see Eqs.(54,92)], /('-dependence 
of the functions (a) Bl=i (even sector), (b) Bl=o (solid line), 
Bl^2 (dashed line) (odd sector), for the trimers. To avoid 
diverging functions, these functions were multiplied by K in 
(a) and by in (b). The normalization is arbitrary. 



sociated by the interaction with the incoming atom and 
the scattering is elastic. Furthermore, a multipolar ex- 
pansion can be performed in terms of the atom-dimer 
relative orbital momentum. In the mathematical limit 
of a vanishing kinetic energy, the atom-dimer incoming 
wave is a s-wave and the scattering is characterized by 
the atom-dimer scattering length Cad that we shall calcu- 
late. To next order of the multipolar expansion the atom- 
dimer incoming wave is a p-wave and we shall calculate 
a corresponding atom-dimer scattering volume V^*^. 

The property of elastic scattering at < rigorously 
holds for the model Hamiltonian Eq.(15), since we have 
shown that it admits at most one dimer state (with rota- 
tional degeneracy). Reality with atoms goes beyond this 
model Hamiltonian: there exist in general deeply bound 
dimers, which can make the atom-dimer scattering in- 
elastic even at arbitrarily low relative kinetic energy. The 
corresponding three-body loss rate is estimated in sub- 
section VB. 

S'-wave atom-dimer scattering: Since we have assumed 
a negative total energy E < Q, there cannot be a free 
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incoming three-atom state so that = in Eq.(38). 
In the center of mass frame, the incoming state is an 
atom impinging on a dimcr with vanishing kinetic en- 
ergy; in the subspace with one atom and one molecule in 
the closed channel, this corresponds to a relative orbital 
angular momentum L = 0, that is to a total momentum 
J = 1 since the molecule is of spin unity. According to 
subsection IV B, the incoming state is in the odd sector. 
Mathematically, this scattering experiment corresponds 
to the following splitting for (3, into the sum of an incom- 
ing wave of zero momentum (a delta distribution in K 
space) and a scattered wave which is a regular function 
of X: 



of the incoming wave for /3, that is in the subspace of one 

1/2 

atom and one molecule, includes a factor p^i^^^^- 

In practice, one injects the form Eq.(44) into Eq.(40). 
The i5(K) term gives a zero contribution in the first term 
of the left hand side, since K = and E = — i^dim leads 
to fcioi = iqdim and l//(fcrei) = 0. The 6{k) inserted in 
the second term of the left hand side of Eq.(40), that is 
the integral term, produces a smooth source term in the 
left hand side. 



T(K)^4vrpJ/l,(e,.K)K 



-5b^K^/8 



(45) 



/3(K)^(2^)35(K)p;;/l,e,+/3°"*(K 



(44) 



This is of the form Eq.(41), the delta being in the L = 
sector. Note that the incoming dimer has a probability 
Pcioscd to be in the closed channel, so that the amplitude 



One is left with a linear and inhomogencous system for 
the vectorial function /3°"' (K) , which is then taken of the 
form Eq.(41), with coefficients Bl%{K) and Bl%{K). 
The explicit form of the resulting system is derived in the 
appendix A, and we obtain 



J 



DiK) 



Bl%{K) 
Bl%{K) 



+ 00 



Bl%{k) 



1/2 
closed 



Idim 



(46) 



r 



where we have introduced the diagonal part 



DiK) 



3/(fc,.el) 



(47) 



and the two by two matrix M{K, k) is given in the ap- 
pendix. 

Let us start with an intuitive presentation of the re- 
sults. We expect that, at low K, the scattered wave in 
the L = channel diverges as so that we set 



(48) 



In position space this indeed corresponds to the large 
r behavior 1 — ftad/?", where r is the distance between 
the molecule and the atom, so that Oad is indeed the 
atom-dimer scattering length. In the channel L = 2, 
the outgoing wave is expected to scale as 1/r^ at large 
r, because of the centrifugal barrier; this corresponds to 
Bl=2{K) having a finite limit in K = 0. 

What typical values of Oad can we expect ? For Vs > 
and much larger than , the scattering amplitude of two 
atoms has a modulus < l/a^cs, which is a small value at 
most of the order of b. For k ~ qdim <C a^cs, one finds 
that |/(fc)| — 1/arcs- One may then expect intuitively 
that Oad weakly depends on Vs, and is at most of the 
order of 1/ares, that is at most ~ b. This expectation 
is correct, see Fig. 5, except close to the threshold for 
the existence of a trimcr in the odd sector, where Cad 
diverges. 



We now turn to a more rigorous analysis of the integral 
equation (46). The key ingredient is the low K behavior 
of the various coefficients for qdim > 0. Consider the diag- 
onal term D{K). As we have already mentioned, D{K) 
vanishes in K = Q: since here k-^ci = *('?dim + 3i4r^/4)^/^, 
we see that fej-ci is an expandable function of K which 
varies to second order in K. The same conclusion holds 
for D{K), which therefore vanishes quadratically in K = 
0; in the limit qdimb ^ 1 we find the simple result 



lim 



DiK) 



(49) 



Consider next the coefficients of the matrix M(i^, fc). 
From the explicit expressions given in the appendix A, 
we obtain for a fixed k: 



lim MiK, k) 

K~*0 



kye 



fc2 



Idim 



1 -1 





(50) 



Assuming that the functions k^Bf^^ik) arc bounded, we 
find that DiK)Bl%iK) has a finite limit in K ^ 0, 
obeying the exact relation 



]im^DiK)Bl%iK) 



— / dk- 



2 

3^ 



^dim 



\B 



L=0 



Bl%]ik). 



(51) 



On the contrary, we find that the second line of the ma- 
trix MiK, k) vanishes quadratically for K ^ 0, and the 
source term also vanishes quadratically in K, so that 
B'^^i^) has indeed a finite limit in K = 0. 



11 



The existence of a well defined limit for Cad in the large 
scattering volume can also be argued in simple terms. All 
the coefficients in the integral equation (46) have a well 
defined limit for qdim 0. In particular, the diagonal 
term in this limit assumes the simple form 



lim D{K) 



4 



h{K)] 



(52) 



where /i(A') = qe'^''' erfc((?6), with q = VSK/2, varies 
monotonically from zero to l/{n^^^b); since arcs > 
l/(7r^/^&), the expression in between square brackets 
cannot vanish. Taking as new functions Gq{K) = 
K^Bl%{K) and G2{K) = K^Bl%{K), one faces for 
Vs +O0 an integral equation of the form 



G(A') - 0[G] = S(A') 



(53) 



where the source term is the infinite limit of the right 
hand side of Eq.(46), and O is a bounded operator de- 
pending on b but not on a^^s- The value Go{0) is finite 
for Vs = +00, so is the atom-dimer scattering length. 
As shown in the appendix, a simple transformation can 
make the operator O hermitian; numerically, one finds 
that the positive part of the spectrum of O consists of 
a continuum extending from zero to l/(47r^/^fe), and of 
a discrete state of energy above the continuum. We see 
that Q!res/4 cannot match an eigenvalue of the continuum, 
but can indeed match the discrete eigenvalue, for 



„odd 



0.69208/&. 



(54) 



This particular value of a^cs corresponds to the thresh- 
old for the formation of an odd trimer at = 00, and 
the corresponding eigenvector was plotted in Fig. 4a. For 
ttics close to the threshold value, the solution of (53) ac- 
quires a diverging component on this eigenvector; since 
the eigenvector has a value Go(0) 7^ in ii' = 0, this 
leads to an atom-dimer scattering length Oad diverging 



as l/(a 

r 



odd 
^th 



The writing (53) also makes it clear that asymptotic 
expressions can be obtained in the narrow resonance 
limit ares^ 2> 1: in this limit, the term proportional 
to Q!res dominates over the bounded operator O, which 
can thus be treated as a perturbation. To leading or- 
der, {Go{K),G2{K)) = 4S(A')/aros, which, injected into 
(51), gives the asymptotic equivalent 



Oad 



32 



(55) 



valid in the limit of large Vg/b^ and large a^^sb- We have 
checked that this relation is obeyed by the numerical re- 
sults. It is important physically to point out that, as we 
shall see in section VI, this asymptotic result no longer 
holds in presence of direct interaction between atoms in 
the open channel. Anyway, it clearly shows that Oad de- 
pends not only on the effective range parameter arcs but 
also on the range 6, which is sensitive to the microscopic 
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FIG. 5: (Color online) Atom-dimer scattering length Oad as 
a function of Ores for fixed values of the scattering volume 
Vs, Vs = W^b^ (black solid line), V, = 10^6^ (red dashed 
line), Vs ~ Wb'^ (green dotted-dashed line). The divergence 
of Oaci coincides with the threshold of existence of a trimer in 
the odd sector. In the limit of a broad Feshbach resonance 
otieab I/tt^^^, Oad tends to 0.2b. 



details of the model interaction. In this sense, the large 
scattering volume limit of Oad is not a 'universal' quan- 
tity. 

This differs from the bosonic case on a narrow Fesh- 
bach resonance, where the atom-dimer scattering length 
is a function of the scattering length a and the effective 
range rg only, as soon as a greatly exceeds the range of 
the potential; furthermore, this function is not bounded 
in the large a limit, but rather exhibits, on top of an over- 
all linear growth with a, a series of divergences for values 
of a/re corresponding to a threshold for the formation of 
an Efimov trimer [58]. 

P-wave atom-dimer scattering: We now assume that the 
incoming atom-dimer relative wave is a p-wave, that is 
it has a unit orbital momentum L = 1. In the subspace 
with one atom and one closed-channel molecule, the cor- 
responding orbital wavefunction is obtained in momen- 
tum space from the Iow-ATq expansion of the Dirac dis- 
tribution corresponding to a molecule of wavevector Kq 
impinging on an atom of wavevector — Kq: 

{2TTfS{K-Ko) = (2^)3 [6{K) - Kq • (grad<5)(K) + ...], 

(56) 

and one may take exactly E = — -Edim at this order. Since 
the molecule has a spin Snio\ = 1 this may correspond to 
a total spin J = 0, 1 or 2. The present work is restricted 
to a total spin J = 1 , and the corresponding ansatz turns 
out to be in the even sector: 

/3(K) = i2nf{-Ko)pl(L, x 

X {[e, • (grad5)(K)]e.y - [e^ • (grad <5)(K)]e J 



K 



(K-e^)e.y - (K-ey)e^]. 



(57) 



Wc insert this ansatz in the integral equation (40), keep- 
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ing in mind that here = 0- The part of the ansatz in- 
volving the gradient of the Dirac distribution gives a van- 
ishing contribution in the diagonal term of the equation 
(since l//(fci-ci) vanishes quadratically in K = for the 



total energy E = — E'dim), but gives a non-zero, smooth 
contribution in the integral term, serving as a source term 
for the scattered wave -B£'^i- Performing the angular av- 
erage as detailed in the appendix A we obtain 



J 



2 /■+°° 

D{K)BZ'tl,iK) ~- dk k^[CoiK, k) - C2{K, k)]Blti{ky 



l^K- \ 1/2 SZ^^K) 



2 ' 
Idim 



(58) 



where the functions Cq and C2 are defined in the ap- 
pendix taking q = (/dim, and D{K) is given by (47). 

The analysis performed for the atom-dimer s-wave 
scattering is readily extended to the present p-wave scat- 
tering. Since the inhomogeneous term in the right hand 
side of (58) vanishes linearly in if = and the diagonal 
part D{K) vanishes quadratically, B^-y{K) diverges as 
\/K. Such a low- iiT behavior was expected: From Eq.(l) 
expanded to first order in the incoming wavcvector, here 
called Kq rather than k, one obtains for the wavefunction 
at large distances 



V-Ko (r) - «Ko • r 



3V 



ad 



(59) 



Taking the Fourier transform with respect to the relative 
atom-molecule coordinates r = rmoi ^ Tat leads to the 
low- if behavior 



^Ko(K) ~ -(27r)3(Ko • grad)5(K) - 12^V, 



ad 



K K 







So we conclude that 



K ■ 



(60) 



(61) 



We thus obtain the exact relation: 



1™ 



ad _ ^ / dk 



^out 



l(fc) 



^dim 



k^ + ql^ 47r(-ifo)p^/ld 
(63) 



Next we take the limit of an infinite scattering volume 
in (58), that is we take gdim — > 0. The source term now 
diverges as 1/K in if = 0; since D{K) vanishes as K^, 
we expect that the function 



Soo(if) 



hm BrMK) 

g— i--|-00 



(64) 



diverges as 1/K^ in if = 0. To check the existence 
of Boo as a limit, one thus has to check that the inte- 
gral in (58) does not have a divergence in = for 
such a 1/k^ behavior of the B{k) function: the fac- 
tor k^ of three-dimensional integration and the fact that 
Co (if, k) — C2{K, k) vanishes linearly with k indeed bring 
an overall k^ factor that compensates the divergence. As 
a consequence it is reasonable to assume that there exists 
a constant C such that 



\Bl%{K)\ < 



(65) 



From the numerical solution of (58), the atom-dimer scat- 
tering volume Vg'^ seems to scale as the atom-atom scat- 
tering volume Vs itself close to the Feshbach resonance. 
So we plot in Fig. 6 the ratio Vg'^/Vg as a function of arcs 
for increasing values of Vs. Another interesting feature 
is the divergence of VJ"^ at the threshold for a trimer 
formation in the even sector. 

The same analytical techniques as in the case of s- 
wave atom-dimer scattering may be used to predict the 
scaling of V|"i with V^. First we divide (58) by if and 
we take the limit if — > 0. As discussed in the s-wave 
atom-dimer scattering case, D{K)/K^ has a finite limit, 
so does D{K)Bf^^^{K) / K . Furthermore one can show 
from (A9) (with q = gdim) that 



lim 



Co(if,fc)-C2(if,fc) 2 



if 



^dim 



(62) 



uniformly in if and Vs- This allows to show that the 
integral term in (63) is 0{l/qd\m) and is thus negligible as 
compared to 1/Qdim- Using (49) and gdim ^ l/(arGsVs)^/^ 
we obtain 



V: 



ad 



:V., 



(66) 



This result corresponds to the dotted line in Fig. 6. 
Strictly speaking, it asymptotically holds for all values 
of a-ccab except right on the threshold for the even trimer 
formation, for reasons that are explained in subsection 
IV E. Away from this threshold we thus reach the im- 
portant conclusion that, very close to the Feshbach reso- 
nance, the atom-dimer scattering volume for a total an- 
gular momentum J = 1 is a 'universal' quantity in the 
sense that it does not depend on the range b of the inter- 
action, but only on the atom-atom scattering volume. 
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FIG. 6: (Color online) Atom-dimer scattering volume VJ"^ for 
a total spin J = 1 (see text) as a function of Qrosb, for a 
fixed value of the atom-atom scattering volume Vs/b^ = 100 
(dashed-dotted green line), Vs/b^ = 10^ (dashed red line), 
Vs/h^ = 10* (sohd black line). To reveal the scaling of V,^"* 
with Va close to the Feshbach resonance, V|''* is expressed in 
units of Vs . Dotted horizontal hne: analytical prediction (66) 
in the limit Vs — > -|-c». 



E. Scattering of three atoms: recombination rate 

In this subsection, we consider the case of three in- 
coming atoms, in the form of plane waves of wavevectors 
k]*, and kg. Without loss of generality wc move to 
the center of mass frame and assume a vanishing total 
momentum. We consider the case of a positive scatter- 
ing volume Vs > 0, so that there exists a dimer state 
in the two-body problem, that may be populated by the 
collision of three atoms. The goal here is to determine 
the rate with which such a dimer state is formed, the 
so-called recombination rate. 

For this physical situation, the total energy is positive 
so Aq in Eq.(38) does not vanish, but rather defines the 
state of the three incoming fermions, 



l*o) 



(fikcPK 



^o(K,k)alj^ 



+k'^iK-k" 



«^-kIO). (67) 



Setting Ko = k? + k§ and ko = (k? - k^)/2, one has 

Ao(K, k) = {2^f5{K - Ko)i [<5(k - ko) - <5(k + ko)] . 

(68) 

To derive a simplified expression in the low incoming ki- 
netic energy limit. 







(69) 



which implies E <C i?dim, one expands Aq in powers of ko 
and A'o. The expression between square brackets gives 
fco times a gradient of delta -l-O(fco). The expansion in 
powers of Ko gives 

A(K,k) = (2^)6[5(K)-Ko-(grad5)(K)+0(A'2)] 
X [-ko-(gradJ)(k)+0(fc3)] . (70) 



One has to keep the leading order in ko ~ Ko giving a 
non-zero value for the incoming state |^'o)- Keeping the 
first term in the first factor gives a vanishing contribution 
so that one has to keep the second term: 

Ao(K,k) ~ (27r)6 [Ko • (grad^)(K)] [ko • (grad<5)(k)] . 

. 

This choice for Ao corresponds to the limit of a vanishing 
total energy, so that we now take E = 0. 

This expression for Ao, when inserted in Eq.(40), gives 
in the right-hand side the source term 



mA 

{[Ko-(grad<5) (K)] ko 



[ko-(grad5) (K)] Kof72) 



Since the gradient of S can be seen as the product of K 
with an isotropic distribution, one finds that this source 
term is in the even sector, of the form Eq.(42), where 
is taken along the direction of Ko and ey is taken along 
the direction of kjj-, the component of ko in the plane 
orthogonal to Ko- We take for /3 the even ansatz with a 
specific form adapted to the present physical situation. 



/3(K)- 

g{K) {[Ko-{gr ad S) (K)] ko 
4TTg{K) 1 



K' - zO+ K 



[ko-(grad5) (K)] Ko} 
[(K • e,)e, - (K • e,)e,] . (73) 



The first term in the right-hand side is motivated by the 
fact that the source term contains a gradient of 6, so that 
f3 has also to contain a gradient of 5. In the second term, 
we have pulled out explicitly a singularity with a pole at 
K = Kdim + *0+, where 



(74) 



which is the value of K given by Eq.(39) when fcioi — 
iQdism keeping in mind that the total energy is here 
~ 0. Physically A'dim is the value i^out of K corre- 
sponding to the motion in opposite directions of a fly- 
ing atom and a flying dimer formed by the three-atom 
collision, and the term iO^ in the denominator of the 
ansatz ensures that this relative motion is a purely out- 
going wave. The conservation of energy indeed imposes 
E = K^^^/4m — Edim, that is Kont — Kdim since we 
assumed E <C i?dim- 

We now inject the ansatz Eq.(73) in the integral equa- 
tion Eq.(40). The bit in gradient of d in the ansatz, when 
injected in the diagonal term of (40), gives a contribution 
which is a distribution of the same structure as the source 
term (72) created by Ao; the function g{K) is adjusted 
to have an exact cancellation: 



g{K) = -(2^)= 



enh^ 3/(fc„ 



771 



A e 



(75) 



where k^^i is defined in Eq.(39) and is equal here to 
i{V3/2)K . When injected in the integral on the left hand 
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side of Eq.(40), the bit in gradient of (5 in the ansatz gives pendix A. one finds an inhomogeneous integral equation 
rise to a smooth function of K (not a distribution). After for g{K): 
lengthy calculations and angular averages detailed in ap- 



- / dK' K'^ [Co (if, K') - C2{K, K')] 
7^ Jo 



9{K') 



-5b^{K^+K'^)/8 



— Aofc^V,- 

niA 



K 



(76) 



where Co and C2 are given by (AlO) and (A12) with k = 
K' and q = 0. It remains to solve this integral equation; 
one notes that there is no delta distribution arising in the 
first term of this equation, since l//(fci.ci) = l//(*'7dim) = 
for K = A'dim, so that iO+ may be omitted in the 
denominator of this first term [59] . 

To obtain the recombination rate from the solution 
g{K) of Eq.(76), we proceed in two steps. First, we cal- 
culate the rate of dimer formation, that is the recombi- 
nation rate, in terms of g{K^i-[a), after having enclosed 
the three atoms in a fictitious cubic box of size L. Sec- 
ond, we construct an operator O such that its expectation 
value in the unperturbed incoming state (67) of the three 
atoms gives the recombination rate; calculating the ex- 
pectation value of this operator for a Fermi sea in the 
thermodynamic limit then gives the recombination rate 
for a macroscopic gas. 

Recombination rate for three atoms: Enclosing the three 
atoms in a arbitrarily large cubic box with periodic 
boundary conditions introduces the following normaliza- 
tion factor in the state vector, 



= of a large atom-molecule separation: 



box 



(rmoi; Tat) ^ ■-Q/9 A e^; 

out L^/'^ r 



(78) 



with r = 7' = r^wi ^ ^a,t is the position of the rela- 
tive particle. The associated probability current for the 
relative particle of reduced mass 2m/3 is then 

j™t = 5~7^^"^ [*;,out5r*7,out] (79) 

7 ' 

^ ^^^\9{K^nn)\' [{er ■ + (e. • e,)^] e.. 

One then calculates the total flux of the current through 
47r steradian and one integrates over the center of mass 
position. Since the flying dimer has a probability ampli- 



tude p 



1/2 

closed 



to be in the form of a molecule in the closed 



channel, it remains to divide the total flux by Pciosod to 
get the rate of dimer formation for three atoms in the 
box, 



I*) 



(77) 



dt 



l5(A'dim)l' 1 



inclosed 



L6' 



(80) 



This is directly seen on the incoming state vector (67): 
each of the three atoms is in a plane wave, with a wave- 
function in the box differing from the free space one by 
the normalization factor l/L"^^^. 

To calculate the probability flux of dimer formation, 
the most convenient is to perform the reasoning in the 
subspace of (35) with one atom and one closed-channel 
molecule, where the formation of a dimer manifests itself 
by an outgoing wave of the molecule of momentum i^dim 
and an outgoing wave of the atom with the same mo- 
mentum in the opposite direction. In momentum space, 
this outgoing wave results from the existence of a pole of 
/3(K) in K = ii'dim, as was made apparent in the ansatz 
(73). In position space, taking the Fourier transform of 
/3(K) and writing g{K) = g{Kdim) + [g{K) - <?(Adh„)], 
we isolate the outgoing wave, and we obtain in the limit 



Recombination rate for a macroscopic gas: To extend 
(80) to a macroscopic number of atoms, we heuristically 
generalize to fermions an operatorial expression derived 
in [60] for bosons: In the bosonic case, the recombina- 
tion rate in a macroscopic gas is expressed in terms of 

-1 3 r ^ -1 3 

■(/'^(R) V'(R) )o, where is the bosonic field op- 
erator and the expectation value (...)o is taken in a 
mean-field state for the bosons not including the short 
range microscopic correlations induced by the interac- 
tion potential [61]. Such a local formula results from the 
assumption that the size of a produced dimer is much 
smaller than the macroscopic correlation lengths of the 
gas, such as the healing length and the thermal de Broglie 
wavelength [60]. 

In the case of fermions, one has to rederive the formula 
since tp^ = 0. This is done in the appendix B and leads 



15 



to the following prescription for the recombination rate: 



dt 



{a,0)e{{x,y).ix,z).iv,^)} 

t / ^\t 



x{{dRj) (OrJ) (^OrJ) (^dR^,4^))o, (81) 

where the field operator and its derivatives are all eval- 
uated in R. This expression involves as a factor the re- 
combination constant /Crcc- not to be confused with the 
recombination rate. In the considered limit of a fcrmionic 
kinetic energy smaller than the dimer binding energy (69) 
we indeed expect /Cj-cc to be a constant, that is not to de- 
pend on the fcrmionic kinetic energy. On the contrary, 
the recombination rate dNtn^/dt will involve a factor 
proportional to the square of the kinetic energy of the 
fcrmions, as predicted in [42] with a different approach. 

To illustrate this point, let us consider the case of a 
spatially homogeneous weakly interacting zero tempera- 
ture Fermi gas. The condition of low kinetic energy is 
then that the Fermi energy ?i^fc|,/(2m) is smaller than 
the dimer binding energy £^dim- It remains to calculate 
the expectation value (. . .)o of Eq.(81) in the Fermi sea 
of the ideal Fermi gas of density n = fc|,/(67r^), to get 



—I 

dt^ 
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(82) 



The factor kp reveals the expected kinetic energy depen- 
dence of the recombination rate. The recombination rate 
will weakly depend on temperature as long as the gas 
remains strongly degenerate, fc^T ^ h^kp/{2m). 
Value of the reconibination constant: We obtain /Croc by 
applying Eq.(81) to our solution of the three-body prob- 
lem. In this case, the uncorrelated state |4'o) over which 
to average in the expectation value (. . .)o is a Slater de- 
terminant with three atoms in plane waves of wavevcc- 
tors k? = ko -I- Ko/2, k^ = -ko -I- Ko/2 and k§ = -Kq, 
respectively. Using Wick's theorem [62], we obtain 

(ko A Ko)2 



dt '^'"^ 



L6 



Equating this expression to Eq.(80) wc obtain 



47r hK, 



dim 



l5(ifdim)P 



9 



(Ko X ko)2pclosed' 



(83) 



(84) 



We find that, as expected, this recombination constant 
does not depend on the incoming energy, that is on 
the norms Kq and ko, in the present limit of vanish- 
ing incoming energy: g{K) is indeed proportional to 
Kok^ = ||Ko X koll, as the source term in the linear 
equation (76) is. 

One solves the integral equation (76) numerically, to 
access g{Kdim)- The corresponding values of the recom- 
bination constant are given as functions of Ofics in Fig. 7a, 
for three values of the scattering volume. As expected, 
a rapid rise of the recombination constant is observed 
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FIG. 7: (Color online) Recombination constant /Crcc appear- 
ing in the expression (81) giving the rate of formation of 
weakly bound dimers when three low energy atoms are col- 
liding, as a function of arcs for a fixed value of the scattering 
volume, (a) /Crcc in units of fdfi /m. The scattering volume 
is, from bottom to top, Vs = lOOb^ (green line), = lOOOfe^ 
(red line) and Vs ~ 10*6'' (black line). Solid lines: numerical 
solution. Dashed lines: asymptotic formula (90). (b) Ratio 
of /Crcc to the asymptotic formula (90), for Vs = 10*6'^. Solid 
line: numerical solution. Dashed line: analytically predicted 
Fano profile (93). The insert is exactly the same figure but 
with a log scale on the vertical axis. 



when one gets closer to the Feshbach resonance, that is 
for increasing values of Vg . For a fixed Vg , one observes 
a smooth dependence of /Crcc with arcs, except in the 
vicinity of Ofres^ = 0.8: both a dip and a peak in /Crcc 
are observed; this singular structure becomes extremely 
narrow in the large Vs limit, both the distance between 
the dip and the peak, and the width of the peak, appar- 
ently tending to zero. These features can be obtained 
analytically as follows, by investigating the large Vs limit 
of (76). 

Let us examine first the diagonal term in the left hand 
side of the equation (76). At low values of K, much 
below 1/6, one can approximate the inverse scattering 
amplitude as 1/ f{iq) ~ -l/{q^Vs)+a ~ a{q^ - qjiJ / , 
where we used Eq.(lO) since wc are close to resonance. 
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Setting q = ^/3K/2, we then see that this diagonal term 
at low energy is close to ag{K)/4:, so it is very smooth. 
At high values of K, of the order of 1/b or larger, one 
can directly set = +00, and one sees that the factor 
of g{K) in this diagonal term decreases smoothly from 
Q;rcs/4 to the positive quantity [arcs ^ l/('''^^^^)]/4 when 
K increases to infinity. 

Let us now turn to the integral term. The value of K' is 
cut to values at most of the order of 1/b by the Gaussian 
factor. For low values of K, below 1/b, an approximate 
expression of the kernel can be obtained [63] , 



Co{K,K')^C2{K,K') 



2KK' 



3(if 2 + if'2) • 



(85) 



This shows that the kernel of the integral part is smooth 
and bounded, even in the low K and K' limit. Neglect- 
ing -R'dim in the denominator of the integral term gives 
a diverging factor 1/K'^ which is however exactly com- 
pensated by the K'^ Jacobian term of three-dimensional 
integration. 

The only source of singularity in the solution g{K) may 
thus be the source term, in the right hand side of (76). 
The presence of a factor Vg will cause g{K) to diverge 
at high Vg; by linearity of the equation, and the 1/K 
divergence of the source will lead to a singular behavior 
of g{K) in if = 0. These two problems can be eliminated 
by taking as unknown function 



F{K) 



Kg{K) 



(86) 



We multiply (76) by K/Vs- The kernel of the in- 
tegral term for F{K) taken in the limit i^dim = 
(thus neglecting K^^^ in the denominator) now behaves 
at low momenta as {K / K')[Cq{K,K) - C2{K,K')] ~ 
(2i^r2/3) / [K^ + K''^) , which remains a bounded quantity. 
Having eliminated the singularity in the source term, we 
can suppose that 



lim F[K) 



(87) 



is a regular and bounded function, |Foo(iir)| < |i^|max- 
Then, if one uses the approximate expression (85), one 
easily sees that the integral term in the equation for 
Fao{K) is bounded by 



2 [+°° , , 2A'2 

- / dR'— — — 

7: J, 3(^2 + 



F\ 



:K\F\. 



so that it tends to zero in if = 0. Since the source term 
for i^oo is non-zero for if = we conclude that 



i^oo(O) 



1447rr 



(89) 



Using the useful expression (33) for the probability to find 
the dimer in the closed channel, we thus obtain analyti- 
cally the asymptotic value of the recombination constant 



for large V,/b^ [64]: 

/Crcc~/C™* = -(48^)2 



3ar 



(90) 



We first note that this result, contrarily to the atom- 
dimer scattering length, is 'universal', that is it does not 
depend on the potential range b but only on the param- 
eters Vs and a entering in the low-fc expansion of the 
two-body scattering amplitude. In particular, (90) is not 
specific to our choice of a Gaussian cut-off function in 
X(k), as we have checked for a general cut-off function 
that is derivable with respect to k. Second, the exponent 
governing the dependence in is the same for a broad or 
a narrow Feshbach resonance. It may thus make sense to 
compare this prediction to the earlier work of [42] , where 
a numerical calculation was performed for a resonant in- 
teraction in a single channel model: the recombination 
rate was found to increase as a power law in Vs, with 
an exponent argued to be equal to 8/3. Since 8/3 and 
5/2 differ by about 6% only, it seems difficult to see this 
difference from the numerical results of [42] . 

To sec how our numerical results approach the asymp- 
totic prediction (90), wc have plotted in dashed lines in 
Fig. 7a the asymptotic value IC^^^™^^, as a function of 
arcs, for the considered values of the scattering volume. 
For increasing values of Vs, we indeed observe conver- 
gence of the ratio /Crcc/A^rcc™''* to unity, but this con- 
vergence is not uniform in aros^ the singular structure 
already apparent in Fig. 7a becomes narrower and nar- 
rower for increasing Vs, but e.g. the peak in this singular 
structure leads to increasing deviation from unity of the 

ratio jCrec/IC^r"'- 

The existence of this singular structure and the depen- 
dence of the recombination rate on arcs within this struc- 
ture can be obtained analytically as follows. First, we 
formally write the integral equation obtained for Foo{K) 
in the limit Vs +00 for a fixed value of arcs^ 



Qfr, 



'-F^iK)^Io[F^]iK) = S{K) 



(91) 



where the source term S{K) is obtained by multiplication 
of the right-hand side of (76) by K/Vs, and Iq is a linear 
operator, given explicitly in the appendix C. Wc find nu- 
merically that the spectrum of Iq consists of a continuum 
extending from to l/{4:7:^^^b), and of one discrete eigen- 
value above the continuum. If one remembers that, from 
(29), arcs > 1/iT^^^^b), it becomes clear that the homo- 
geneous equation obtained by replacing S with zero will 
admit a non-zero solution uo{K) only for arcs = ^t™'^, 
where mathematically a^^^'^/A is the discrete eigenvalue 
of Iq, and physically a^J^^''" is the threshold value of ares 
for the existence on resonance of an even trimer. The op- 
erator appearing in (91) is indeed the infinite scattering 
volume limit and the zero energy limit of the operator 
L{E) of the subsection IV C on trimers, restricted to the 
even sector. Numerically we find 



„ even 



0.81408/6. 



(92) 
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In presence of the source term S, and for a value of arcs 
slightly deviating from a^™'^, one realizes that a compo- 
nent of Foo{K), proportional to uq{K), may diverge as 



i/K 



^th 



The appearance of such a small denom- 



inator implies that Foo{K) is not uniformly boimdcd in 
arcs in the vicinity of 0;°™", so that the asymptotic law 
(90) may not hold uniformly in arcs- 

This very simply reveals that the singular structure 
in the recombination coefhcient is a consequence of the 
existence of a weakly bound trimer. Quantitatively, as 
shown in the appendix C, by going beyond the Vs = 00 
approximation, one can calculate analytically the con- 
tribution to F{K) which becomes large for arcs close to 
Q,even_ 'jj-^jg ^gg^dg to a Fauo profile [65] 



/Crc 



t^asympt 



(arcs - "o)^ 



(93) 



,^rcs ai)2 4-Aa2 

with the following low-gdim expansions, 

ao air-Z.2ql^h (94) 

(95) 
(96) 



^th - 
ai ^ a?r + 6.294 g2.^6 

Aa ~ 35.89 gj^irrr^^ 
Aa 16 gdim 



ai - ao 



3^3 a\ 



even 
th 



(97) 



where (/dim — 1/ \Jof^^°^Vs- This is in agreement with the 
numerical results at finite but large V^, see Fig. 7b. 

Note that such a Fano profile in the recombination 
constant as a function of the width of the Feshbach res- 
onance does not occur in the case of bosons in the large 
scattering length limit: In the bosonic case, when the 
scattering length a becomes much larger than the range 
and the effective range of the interaction potential, the 
recombination constant, apart from an asymptotic fac- 
tor, has only a bounded oscillatory behavior as a function 
of a or of the potential range [66, 67]. This is a conse- 
quence of the fact that trimers of bosons always exist for 
large enough a, whatever the width of the Feshbach res- 
onance, and they exist in arbitrarily large numbers for a 
arbitrarily large, the so-called Efimov effect; the oscilla- 
tory behavior in the recombination constant then results 
from the successive entrances of new Efimov trimers, as 
a grows. 

Finally, turning back to the fermionic case, we note 
that, in reality, there may exist deeply bound dimers 
even in the limit |Vs| 3> h^, which may be formed 
by the collision of three incoming atoms, in competi- 
tion with the weakly bound dimer. This effect, beyond 
our model Hamiltonian, is discussed in subsection VC, 
where the corresponding recombination constant towards 
deeply bound dimers is estimated. 



support only a weakly bound dimer in the two-body prob- 
lem, that is with an energy much smaller than /mb^. 

In real experiments, with alkali atoms, the van der 
Waals interaction potential is very deep and support 
several deeply bound dimers in the two-body problem. 
These deeply bound dimers can be formed by three-body 
collisions, and liberation of a huge binding energy, lead- 
ing to particle losses. As a consequence, the trimers will 
acquire a finite lifetime, the atom-dimer scattering will 
not be purely elastic, and the scattering of three atoms 
will lead to recombination not only to the weakly bound 
dimer, but also to the deeply bound dimers. 

A first possibility to estimate these inelastic contribu- 
tions is to modify the Hamiltonian, so as to have deeply 
bound dimers in the model, e.g. by including a separa- 
ble interaction potential in the open channel [68, 69], 
or by considering a true potential and introducing the 
adiabatic potential curves in hyperspherical coordinates 
[70, 71]. This however modifies the mathematical struc- 
ture of the problem and is beyond the scope of the present 
paper. 

Fortunately one may easily estimate the loss rate, that 
is the rate of formation of deeply bound dimers, by the 
following recipe, expected to be accurate within an un- 
known approximately constant factor [39, 72]: 



r -^p 

^ loss ~ , 9 <6 



(98) 



where P<;, is essentially the probability that the three 
atoms be all within a volume of the order of b^. More 
precisely, since we are using a two-channel model, this 
probability can be split in two contributions, one coming 
from the purely atomic component (all three particles in 
the open channel). 



P^\= I d3^ldV2d3r3|vI/^'"(rl,r2,r3)|^ 

J p<b 



(99) 



where p is the hyperradius (B2), and the other contri- 
bution coming from the molecular component (with one 
open channel atom and one closed channel molecule). 



|r„ol-rat|<b 



dVmoldVat|*Zr(rmoi;rat)|' 



(100) 

The wavefunctions are here normalized, hence the apex 
"norm" , as we shall explain case by case. It remains 
to apply a Fourier transform to (3 and to an anti- 
symmetrized version of A in (35) to calculate numeri- 
cally the corresponding wavefunctions ^'{Jr™™ and ^f""™. 
But one can also have analytic estimates, by approximat- 
ing the wavefunctions by their small-radius expansions 
[r; — Tj I ^ &, as we shall see. 



V. ESTIMATION OF LOSSES DUE TO DEEPLY 
BOUND DIMERS 



Lifetime of the trimers 



Very close to the Feshbach resonance, the model 
Hamiltonian that we have considered in this work can 



In the case of trimers, the state vector can be nor- 
malized in the center of mass frame, as detailed in the 
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appendix A. One may then calculate the probabilities 
in Eqs. (99,100) numerically; the corresponding rate Fjoss 
then represents the inverse lifetime of the trimer due to 
spontaneous decay into a deeply bound dimer and a free 
atom. We recall that (98) contains an unknown numeri- 
cal factor that depends on the microscopic details of the 
interaction, so the values of the lifetimes that we shall 
obtain are only indicative. 

As a consequence, it seems more interesting physically 
to obtain the scaling laws of the trimer lifetime close to 
the trimer formation threshold, that is when the trimer 
binding energy is ?i^/(m&^): Does the trimer decay 
rate tend to zero on the threshold ? Even if this is the 
case, one cannot immediately conclude that the weakly 
bound trimers are long-lived, because their binding ener- 
gies also tend to zero on the threshold. One rather has 
to see if the decay rate tends to zero faster or not than 
the binding energy E^^^^ of the trimer. To this end, we 
form what we call the quality factor of the trimer: 



Q = 



iTibind 
trim 



loss 



(101) 



This quality factor is shown as a function of qtrim hi Fig. 8 
for an infinite scattering volume, that is in practice for 
l/l^resVsl ^ 9trim- ^^^^ ^^^^ ^^e quality factor Q 

tends to zero at the threshold for trimer formation, which 
is not a positive result. The odd sector is however much 
more favorable (keeping in mind that the quality factor in 
the even sector was multiplied by a factor 20 for clarity 
in the figure): Values of Q much larger than unity are 
obtained already for moderately small values of gtrim^- 
This is due to the fact that Q vanishes more slowly in 
the odd than in the even sector: on the figure, Qodd 
seems to vanish linearly whereas Qeven seems to vanish 
quadratically. 

The scaling of the quality factor with gtrim close to the 
trimer formation threshold can be obtained analytically 
from the low-K dependence of I3{K), considering again 
an infinite scattering volume. 

In the even sector, we have seen in Fig. 4a that the 
function Bl=i{K) right on the threshold diverges at low 
K as 1/K only. This can be shown, as done in appendix 
C, using the fact that the kernel {Co — C2){K, k) vanishes 
linearly in whereas the diagonal part of the integral 
equation vanishes quadratically as for gtrim = 

0. As a consequence, the function B{K) = KBl=i{K) is 
bounded, and the unnormalized state vector \^) of the 
trimer is square integrable in the center of mass frame, 
as is apparent on Eqs.(A25,A26). We thus face the same 
phenomenon as in the two-body case: at threshold, the 
normalized trimer wavefunction is non-zero. As a conse- 
quence, the probability of finding the particles at relative 
distances less than h tends to a non-zero value at thresh- 
old. Floss does not vanish. 



hm Fr™">0, 

qtrimfa— '0 



(102) 



In the odd sector, the situation is more favorable. As 
we have seen in Fig. 4b, for the arbitrary normalization 
chosen in that figure [73], the function Bl=q{K) right 
on threshold diverges as [whereas the Bl=2{K) is 

0{1/K)]. As a consequence, the functions /3(K) and 
thus ^'moi ai'c not square integrable, and (^'atj^at) is 
also infinite, see Eqs.(A27,A28). This means that the 
normalization factor Aft linking the correctly normalized 
state vector to the unnormalized one j^*). 
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vanishes for qtrimb tending to zero. As a consequence, 
P<h and Floss also vanish. 

We now determine the corresponding scaling law. The 
functions -Bl=o and -Bl=2 solve a homogeneous integral 
equation, corresponding to the homogeneous part of (46) 
[that is with the source term set to zero in the right- 
hand side] written for an energy E = —h^q^^-^^/m. One 
can then recycle the reasoning performed below Eq.(46). 
At if <C 1 /6, the diagonal term 



D{K) ~ a,UK^ 
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(104) 



where -R'trim = 2gtrim/V3. On the other hand, the first 
line of the matrix M{K, k) has a non-zero limit for K 
0, see (50). For a choice of normalization of the S's such 
that the functions fc^|i?L=o,2|(fc) are uniformly bounded 
for Ktrim tending to zero, we reach the form 
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where the envelope function £{K) has a finite but non- 
zero limit m. K = Q, and is uniformly bounded as a func- 
tion of K and -ftTtrim [74] . Since the second line of the ma- 
trix M{K, k) vanishes for if ^ 0, we find that -Bl=2 is 
dominated by i?L=o at low K, iiTtrim and can be neglected 
[75]. Inserting the form (105) in the normalization inte- 
grals (A27,A28) we obtain the asymptotic results in the 
gtrim ^ hmit: 
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where p^ioscd ~ 0.185 . . . here at the odd trimer formation 
threshold. This shows that |A/t|^ scales as gtrim^- In cal- 
culating the probability P<b to have the particles 'inside' 
the interaction potential for a non-zero but small gtrim, 
we can take directly the unnormalized wavefunctions ^' 



for gtrim = 0, so that P^b scales as |A/t|^ and 
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and the quality factor Qovcn tends to zero as (gtrim^)^ 



The quality factor Qodd thus vanishes linearly in qtrim^- 
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FIG. 8; For an infinite scattering volume, quality factor of the 
trimers, as a function of gtrimfo, for the even sector (solid line) 
and for the odd sector (dashed line). For clarity, the quality 
factor in the even sector was multiplied by 20. The quality 
factor is defined as the ratio of the trimer binding energy, here 
equal to i?trim = Ti^qtiim/iTT- since Vs = oo, and of h times the 
spontaneous decay rate of the trimer due to the formation 
of deeply bound dimer and a free atom as estimated by the 
simple recipe (98). 



B. Losses in atom-dimer scattering 

We now estimate the loss rate in the collision of an 
atom with a weakly bound dimer. We enclose the atom 
and the dimer in a fictitious cubic box of volume L^, so 
that the state vector in the box can be normalized to 
unity. The box has a size L 3> 6, |aad| so that this state 
vector differs from the free space one by a normalization 
factor only. 
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One has then indeed, in the subspace with one atom and 
one molecule at large distances : 
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and in the subspace with three atoms, for the position 
of the third atom going to infinity for fixed positions of 
atoms one and two [56, 76], 
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so that, after spatial integration of the modulus square of 
these two wavefunctions in the box, using J <Pr |0(r)p = 
1 — Pciosed, one finds that |vl/"°'^"i) is normalized to unity 
in the cubic box, to zeroth order in Qg^^/L and b/L. 

To link this calculation with an experimentally relevant 
quantity, we consider a low density mixture of A^at atoms 
and iVdini dimers in a volume . The loss rate will be 
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The loss constant i^ad is related to Eq.(98) by setting 



A"at = -^dim = 1 in the above equation 



loss 
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which can be checked to be independent of L^. The re- 
sulting atom-dimer loss constant is plotted in Fig. 9 as a 
function of arcs- We see that, for Vs ^ b^, 
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except close to the trimer formation threshold where, 
within the simple recipe, ii'ad diverges. We also see a 
drop of -ftTad to a smaller but non zero value in the limit 
of broad Feshbach resonances, equal to hb/m within a 
numerical factor: this drop is due to the fact that P™^"' 
tends to zero in this limit, so that P<t reduces to the 
atomic contribution -P^'^, which is elsewhere dominated 
by P'^t- 

The property (114) can be understood analytically. 
E.g. to estimate P^^^, one can approximate in 
(100) by its value in Tmoi = fat, which is generically non- 
zero for the odd ansatz. Taking = oo gives a finite 
value for the wavefunction, since i3£'^g(A') diverges as 
l/K^, see (48), and this is integrable on a vicinity of 
K = in three dimensions. This explains the weak Vs de- 
pendence of if ad for large scattering volumes. For arcs^ of 
the order of unity, Eq.(114) then holds from dimensional 
analysis, apart from the divergence close to the trimer 
formation threshold. At large aresb, we see from (53) that 
the scattered wave B£'^o 2i^) is 0(l/ares) and is domi- 
nated by the contribution of the incoming wave oc 6{K.). 
We thus find again (114), with a dominant contribution 
P™j°' from the molecular sector and an atomic sector con- 
tribution which is about 1 ~ p^i^scd ~ 0[1/(q!ics^)] times 
smaller. 

To make the discussion more complete, we also esti- 
mate the rate of formation of deeply bound dimers when 
the atom and the dimer scatter in the p-wave, each with 
a momentum of modulus Kq. In this case the wavefunc- 
tion for distances between the particles <C l/^l^o is ob- 
tained from the ansatz (57) , containing the overall factor 
A'o. On the contrary, if one takes a box size L ^ 1/A'o, 
the normalization factor linking Ivp"™™) to l^") remains 
As a consequence, the p-wave loss constant 
-^■^Tioss will be proportional to Kq. Further- 
more, to estimate P^§\ we can expand the wavefunction 
'S'moi(i"moi; Tat), which is now in the even sector, to lead- 
Tat, that is to first order, which 



Atd = 



ing order in r = rmoi — 
amounts in the Fourier transform of /3(K) to replacing 
exp(iK • r) with iK • r. Taking directly the limit of an in- 
finite scattering volume, we see from (65) that this first 
order estimate remains finite, because K/K^ — 
is integrable around AT = in 3D. As a consequence, 
we expect K^^ to be of the order of TiK^b^ /m, except 
close to the even trimer formation threshold, where it 
diverges within the present formalism. This expectation 
is confirmed by the numerical calculation (not shown). 
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FIG. 9: Loss constant K^,^ due to formation of deeply bound 
dimers in the atom-dimer collision, as estimated by the recipe 
Eq.(98) and the relation Eq.(113), as a function of arcs for 
Vs = 10^6^ (black solid line), Vs = lO^fe^ (red dashed line), 
green dotted-dashed line), ifad is in units of 
hh/m. A divergence of i^ad occurs at the threshold of odd 
trimer formation. 



second term of the ansatz, we take the large scattering 
volume limit, away from the trimer formation thresh- 
old, so that we neglect i^^im ^hc denominator and 
we take g[K) ~ VsFao{K)lK, as discussed around (87). 
This gives again a contribution proportional to Vs- The 
Fourier transform cannot be calculated analytically, but 
we only need the wavefunction for r ^ 6 so that we can re- 
strict to a small-r expansion of the atom-molecule wave- 
function. The linear order in r is the first non-zero one, 
and we obtain the contribution from the molecular sector 
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The expression in between square brackets is not a slowly 
varying function of a-ccs^, because it diverges in the vicin- 
ity of the trimer formation threshold (an artifact of the 
approximations performed here on the second term of the 
ansatz (73)). We have checked numerically that the ex- 
pression can be well approximated by the fitting formula 



Physically, this means that, away from the even trimer 
formation threshold, will be smaller than i^'ad by a 
factor of the order of [Kob)^ ^ that is at least by a factor 
/Vs since we assume here Kq <^ (/dim- 



C. Recombination to deeply bound dimers 

To complete this section, we now evaluate the rate of 
formation of deeply bound dimers in the collision of three 
asymptotically free atoms. We again enclose the atoms 
in a fictitious cubic box of size L with periodic boundary 
conditions, so that the normalized state vector in the box 
is related to the free space one in the large L limit by (77). 
Reproducing the reasoning of the appendix B leading to 
the prescription (81), we can now define a recombina- 
tion constant towards deeply bound dimers, that we call 
IC^^c"^- Applying the resulting prescription to three atoms 
in the box, we obtain the equivalent of (83) for the rate 
of formation of deeply bound dimers, with /Croc replaced 
with /Cfcc^- the other hand, the rate of formation of 
deeply bound dimers is Fioss defined by the recipe (98). 
Equating the two expressions of this rate leads to 
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We consider first the molecular sector. One has to 
perform the Fourier transform of /3(K) in Eq.(73) to ob- 
tain the atom-molecule wavefunction. The Fourier trans- 
form of the first term of the ansatz can be calculated 
exactly: Since it is composed of gradients of delta distri- 
bution in momentum space, it gives a wavefunction vary- 
ing linearly with the coordinates of r = rmoi ~ Tat ■ It is 
found that this wavefunction is proportional to V^. In the 
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The same procedure can be applied in the atomic sec- 
tor. To estimate the purely atomic wavefunction, we ex- 
pand it to leading order in the interatomic distances. The 
zeroth and first order vanish, and we get to second order 

vl/at(ri,r2,r3)c.^ j ^!g!^^(K, k)(x A y) • (k A K) 

(118) 

where we have introduced the Jacobi-like coordinates x = 
ri — Y2 and y = ra — (ri -I- V2)/2. This leads to 
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where the factor /siow depends only on a-^csb, it is of the 
order of 0.005 for a broad Feshbach resonance and it in- 
creases by a factor ~ 5 from broad to narrow Feshbach 
resonances. 

As a consequence, our estimate of the recombination 
constant towards deeply bound dimers, away from the 
trimer formation threshold, scales as 
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for a given arcs, in the large scattering volume limit. We 
thus see from (90) that the formation of weakly bound 
dimers wins over the deeply bound ones in this limit. 
Note that the estimate of /Cre'?^ ^l^o holds on the neg- 
ative scattering volume side of the resonance, still re- 
stricting to the low relative incoming atomic momenta 
A:«l/(arcs|V.|)i/^ 
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VI. EFFECT OF A NON-RESONANT 
INTERACTION IN THE OPEN CHANNEL 

In real life there exists an attractive van der Waals 
interaction between atoms in the open channel, respon- 
sible for a residual interaction in the p wave even very 
far from the Feshbach resonance. This residual interac- 
tion may be characterized by the so-called background 
scattering volume V^*'. Usually it is assumed that this 
residual interaction is weak, |V^^| « 6^ where the in- 
teraction range b is of the order of the van der Waals 
length, so that it is neglected in the vicinity of the Fesh- 
bach resonance as compared to the effect of the coupling 
to the closed channel [30, 36]. However, with the pure 
closed channel coupling Hamiltonian (15) used in this pa- 
per, we found that several quantities were depending not 
only on the low-fc scattering properties parameterized by 
Vs and a, but also on the range b of the potential, such 
as the threshold for trimcr formation, which raises the 
issue of their dependence with the microscopic details of 
the model. Furthermore, we found that the atom-dimer 
scattering length assumes values smaller than b for broad 
Feshbach resonances, so that it is not evident that the 
residual interaction is really negligible. 

To address these questions, we model the residual in- 
teraction by a separable potential of coupling constant go 
with the same cut-off function x(k) as in the closed chan- 
nel coupling. This amounts to adding to the Hamiltonian 
(15) the open channel interaction 
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We determine by relating it to V^^ from the solu- 
tion of the two-body problem. Then we solve the three- 
body problem again with the simultaneous inclusion of 
the closed channel coupling and the open channel inter- 
action. 



A. Modification of the two-body problem 

The calculations proceed along the lines of subsection 
niB. The same ansatz (19) for the two-body state vector 
applies; the new term emerging from the action of V^pon 
is simply 

/d^k 
^X(k)4alkl0) (122) 



where we have set 
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From Schrodinger's equation at energy E, we find the 
remarkable property 



(ii;-i?„oi)^ + 7-o. 



(124) 



so that 7 can be expressed in terms of (3 and the new 
reduced scattering amplitude has a form very similar to 
the previous one (24), 
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where E = h^kQ/m is the energy of the two-body scat- 
tering state. This leads to a modified expression for the 
scattering volume. 
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We first analyze the result very far from the p-wave 
Feshbach resonance. Taking the limit E^^i ^ oo in the 
above expression gives the background scattering volume 
as a function of the open channel coupling constant. 
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Since the van der Waals interaction is attractive, we take 
5o < in all what follows. Then one sees that the back- 
ground scattering volume has a dependence with |(7o| sim- 
ilar to the left part of the Fig.l calculated for a square 
well potential: For increasing values of m\gQ\/h'^ starting 
from zero, V^^ decreases from zero to — oo, it diverges on 
the critical value 
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then it decreases from -|-oo down to 2tt^/%^. The diver- 
gence is due to the formation of a dimer in the open chan- 
nel, and this dimer is deeply bound when g^ is away from 
the critical value ^q. The existence of a deeply bound 
dimer would deeply change the physical nature of the 
three-body problem with respect to our previous analy- 
sis: E.g. it would open a decay channel to the trimer, 
which would not exist as a true stationary state anymore 
but at most as a resonance. We thus take from now on 
l.9o| < 1.90 1 so that V^^ < 0. We shall keep in mind that 
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for a non-resonant interaction in the open channel. To 
be complete, we have also calculated the value of the 
parameter a in the presence of open-channel interactions 
only: 
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which can have any sign since the open-channel interac- 
tion is not resonant. 

We now come back to the vicinity of the Feshbach res- 
onance, where the closed channel coupling is no longer 
negligible. First we can prove that, under the condition 
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V^^ < 0, there is no bound state in the two-body prob- 
lem for v., < and there is one for Vg > 0, as expected; 
the proof was obtained using the argument of a mono- 
tonic variation of an appropriate function as done in the 
paragraph below (30). 

We restrict for simplicity to the exactly resonant case 
Vs = oo. The corresponding values of the closed-channel 
detuning E^oi and of a are then: 
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We see that the range of variation of a-^cs remains the 
same as in our previous model. The scattering ampli- 
tude, analytically continuated to negative energies E = 
—h^q^/m, q > 0, can be put in the simple form 
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(133) 

making it apparent that the main effect of the open chan- 
nel interaction is to introduce a new scale gopcn for the 
wave vectors, such that 
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This allows to reach first conclusions on the effect of the 
open channel interaction on the properties of the original 
model (15): 

• for a broad Feshbach resonance A ^ h'^b^/'^/m, we 
find that qopenb > 1, so the open channel inter- 
action should have a weak effect. In particular, 
this suggests that the trimer states in the regime 
of rather broad resonances should be weakly af- 
fected. Furthermore, if V^^ is weak enough to have 
qopenb > 1 at the threshold for the formation of the 
trimers in our previous model. 
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then the threshold itself should be weakly affected 
by the open channel interaction. For a narrow Fes- 
hbach resonance, one has (jopon^ < 1 (except for 
|V^^| <C b^); in this case the effect of the open 
channel coupling is more difficult to guess: It may 
depend on the considered quantity and a more de- 
tailed analysis is required. 

• if Qopcn was smaller than the estimate gdim ~ 
l/(Q!resVs)^/^ for the wave vector associated to the 
dimer binding energy close to the resonance, then 
the open channel interaction would have a dramatic 
effect. But one finds qli^n/llpen < l^s^/^sl so that 
<Zdim *C gopcn in thc rcsonant regime. 



A last relevant quantity is the probability to find the 
dimer in the closed channel. In the limit Vg +oo, 
after some calculation, we find 
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The relation (33) is therefore affected by the open channel 
interaction, whereas (34) is not. 



B. Modification of the three-body problem 

We now solve the three-body problem in presence of 
both the closed channel coupling and the open channel 
interaction. The previous ansatz (35) applies, but a new 
term arises in Schrodingcr's equation of eigenenergy i?, 
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Schrodingcr's equation projected onto the molecular 
subspace is unaffected by the open channel interaction, 
since there is only one atom in that subspace. So Eq.(36) 
still holds exactly. Then one immediately sees that a 
simple relation relates 7 to /3: 



E-E, 



mol 



Am 



2A 



-7(K) = 0. 



(140) 



One may then take as unknown any convenient combi- 
nation of 7 and /3. 

In Schrodingcr's equation projected onto the atomic 
subspace, a new term appears, but of the same structure 
as the term involving /3 in (37). Thus the modification 
to (38) is minor. 
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This results immediately suggests which combination of 
7 and (3 is convenient: We introduce 

/3,ff(K) = /3(K)-|7(K), (142) 

the overall factor being such that /3gff reduces to f3 in the 
absence of open channel interaction. As a consequence, 
our unknown field is now 
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Eliminating f3 and A in terms of this unknown field in 
(36), we find after some calculations and using (125) that 
/3(,ff (K) solves the same equation (40) as (3 in our previ- 
ous model, provided that the modified scattering ampli- 
tude (125) is used. As a consequence, all the numerical 
and most of the analytical techniques developed for the 
previous model may be reused for the new model. 

Existence of weakly bound trimers: We reproduce the nu- 
merical calculations of subsection IV C with the scatter- 
ing amplitude modified by the open channel interaction. 
We restrict for simplicity to an infinite scattering volume: 
See Fig. 10 giving qtiimb as a function of ares&, where the 
energy of the trimer is —h'^q^^^^/m, for (a) the even sector 
and (b) the odd one. Then we see that the trimer state 
still exists for low a^esb- For a small background scatter- 
ing volume V^^ = —6'^, its energy dependence with ttres^ 
is only weakly affected by the open channel interaction, 
as expected from the qualitative condition (135), and the 
threshold is only slightly shifted. For a much more neg- 
ative background scattering volume Vj's = — lOfe'^, there 
is simply a larger shift in the odd sector, but the conclu- 
sion is radically changed in the even sector, see Fig. 10a: 
the trimer seems to exist now for all values of arcs 6- This 
may be understood as follows: in the large ares^ limit, the 
coupling A to the closed channel tends to zero, so does 
-^moi' (131), and for a fixed non-zero q, the scattering 
amplitude (125) converges to the one of a single channel 
model with a scattering volume V^^. If |V^s|/6^ is large 
enough, then this single channel model can indeed sup- 
port a trimer [77]. On the other hand, the fact that the 
threshold for trimer formation survives in the odd sector 
up to higher values of |Vj8|/6'^ than in the even sector 
may be understood from the qualitative argument (135). 
S'-wave atom-dimer scattering: We consider for Vs > 
the scattering of an atom on a dimer in the limit of a van- 
ishing relative kinetic energy, so that E = — i?dim, where 
Edim is the dimer binding energy for the new model, and 
the corresponding s-wave scattering is characterized by 
the atom-dimer scattering length Oad- Then the field 
/3(K) is given by the ansatz (44). As a consequence, the 
effective field (143) will have the same structure; the delta 
distribution (5(K) will simply be multiplied by the factor 
1 - go (-Edim + £',noi)/(2A2 ). Also the part of /3^fj diverg- 
ing as will be related to the part of /3(K) diverging 
as 1/K^ by exactly the same factor. So that one may 
read the value of the atom-dimer scattering length di- 
rectly from the effective field. As a consequence, one has 
to solve the same integral equation (46), just changing 
the diagonal part D{K) to account for the new scatter- 
ing amplitude. The corresponding numerical results for 
Oad are presented in Fig. 11, as functions of ares, for an 
infinite scattering volume Vs. For the considered values 
of |V]'s/53|^ ^jjg open channel interaction does not qual- 
itatively change the result: There exists a threshold for 
the formation of a trimer in the odd sector, see Fig. 10b, 
and we recover the divergence of a^d at this threshold 
already observed in our first model. 

What happens to the atom-dimer scattering length in 
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FIG. 10: (Color online) In presence of open channel interac- 
tions, characterized by a fixed value of the background scat- 
tering volume Vj^, parameter gtrim of the trimer (when it 
exists) (a) in the even sector and (b) in the odd sector, as 
a function of arcsb. Here the scattering volume Vs is infinite 
and Qrcs is the corresponding value of the parameter a for 
the modified scattering amplitude, see Eq.(132). Solid line 
(black): = 0. Dashed line (red): V,'"^ = -b^. Dashed- 
dotted line (green): V^'''^ = -lOb^. 



the large ares^ limit ? Is the analytical prediction (55) 
obtained in our previous model still valid ? In presence of 
interactions in the open channel, it seems surprising that 
flad can tend to zero at large ares^, since the scattering 
amplitude for A for a finite q tends to the one of 
a single channel model with a scattering volume V^^. 
This limit however is not reached uniformly in q: e.g. 
for q ~ (?dim, q ^ (Zopcn and the scattering amplitude 
remains very close to the one of the two channel model 
with no open channel interaction. This non-uniformity 
of the A — > limit is also revealed at Vs = oo from the 
fact that the parameter arcs is very different from the one 
(130) that one would have in the absence of coupling to 
the closed channel. Mathematically, Eq.(51) still applies, 
if one considers the effective field (3^ff(K) rather than 
/3(K), but in the infinite scattering volume limit analysis, 
the function h{K) in (52) is changed by the open-channel 
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FIG. 11: (Color online) Atom-dimer scattering length Oad as 
a function of Qrcsb for the model including the open channel 
interaction, for an infinite scattering volume Vs. Solid line 



(black): Vs*5 = 0. Dashed line (red): 
dotted line (green): V^*^ = -106^. 
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with q ~ ^/3K/2. Contrarily to our previous model, the 
function h{K) increases from to arcs when K increases 
from to infinity, in practice to values ^1/6. We can no 
longer assume l/(Q;ros — ^{K)) ~ 1/ares for all K in the 
large arcs limit, so we have only the weaker result that 
the right hand side of (51) is 0(1) in this limit. Since 
h{K = 0) = 0, D{K)/K^ still converges to at.os/4 in 
K = and we get 



flad = O 



(145) 



We have successfully compared this analytical predic- 
tion to the numerics. For very large values of aicsb 
(not shown), we numerically find for V^^ < that Oad 
tends to zero as C/aros, where the constant C depends 
on the background scattering volume, C ~ —0.69 for 
V^g = -b^, and C ~ -4.1 for V^'s = -lOb^. 
P-wave atom-dimer scattering: Now the incoming atom 
and dimer have a relative orbital momentum L = 1, 
so the vanishing kinetic energy limit of the scattering 
E —Edim is characterized by an atom-dimer scatter- 
ing volume V^'^, with the ansatz (57) for /3(K) in the 
sector of total spin J = 1 (see subsection IV D). Since 
the factor linking jS^g^K) to /3(K) varies only quadrati- 
cally with K, it may be replaced by its ii' = value in 
front of the gradient of 6(K), so that /3gfj and (3 have 
the same low-i^T behavior, from which VJ"^ is readily ex- 
tracted. We numerically solve (58) updating the values 
of the scattering amplitude f{k) and the dimer binding 
energy i?dim to include the open channel interaction. The 
dependence of VJ"^ with arcs is shown in Fig. 12, for a fixed 
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FIG. 12: (Color online) In the model including an open chan- 
nel interaction, atom-dimer scattering volume V^'^ for a total 
spin J = 1 (as detailed in subsection IV D) as a function 
of Urcsb, for a fixed value of the atom-atom scattering volume 
Vs = W*b^ (solid black line), but for various background scat- 
tering volumes in the open channel: V^^ = (black solid line), 
Vjs = -foS (dashed red line), V^*^ = -lOb^ (dashed-dotted 
green line). VJ'^ is expressed in units of Vs. Dotted horizontal 
line: analytical prediction (66) in the limit Vs +oo. 



and large scattering volume and for various values of 
the background scattering volume. It is apparent that 
is weakly affected by the open channel interaction, 
apart in the vicinity of the even trimer formation thresh- 
old (if it exists). This can be understood analytically, 
realizing that the reasoning leading to (66) still applies 
in presence of the considered open channel interaction. 
This confirms the 'universality' of the asymptotic behav- 



v: 



ad 



Vs 

3 ' 



(146) 



which was expected in subsection IV D from the fact that 
it does not depend on the potential range b. 
Recombination rate to weakly bound dimers: Finally, we 
consider the scattering problem of three atoms in the 
zero total energy limit, for > and in presence of 
open channel interactions. We take for the effective field 
/Sgff (K) the same ansatz as in (73), putting the subscript 
"efF' on the various functions of the ansatz. The func- 
tion Gcs{K) is still given by (75) but the scattering am- 
plitude is now changed. As GcS is however multiplied 
by the gradient of a delta, only its first derivatives in 
K = matter, so only the scattering volume comes out, 
and gos obeys the integral equation (76) with the same 
source term and the updated scattering amplitude. In 
the large scattering volume limit, away from the even 
trimer formation threshold if it exists, we can thus recy- 
cle (89), since it involves only the low- momentum behav- 
ior K < A'dim ^ (Jopcn of the scattering amplitude and 
properties of the function Cq — C2 ; this leads to 



ff eft (A'dim) 



1447r?i^ 
mAarcs Kdu 



(147) 
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On the other hand, the reasoning leading to the recombi- 
nation constant /Croc in (84) still holds, with the function 
g{K) (and not gcs{K)) and the closed channel proba- 
bility Pciosod modified by the open channel interaction. 
The functions g and geff differ, because f3 and f3^ff dif- 
fer by a factor depending on K. Since the factor in 
between square brackets in (143) cannot vanish for our 
choice < 0, f3{K) and (3^^ both have a single pole in 
K = A'dim, with 

ff(ifdim) = gMKdin.) (^43^ 

In the large Vs/b^ limit, one can neglect gdim in the de- 
nominator of this expression, which amounts to neglect- 
ing qdim with respect to gopcn- Then we use the expres- 
sion (136) and remarkably we recover exactly the same 
asymptotic behavior for /Crcc as in the previous model, 
under the assumption V^^ < 0: 

IC..J^>^^'^iASnr(^Y\ (149) 

This indicates some 'universality' of this result, which 
could be hoped from the fact that it does not depend on 
the interaction range b. 



VII. CONCLUSION 

We have solved the free space three-body problem for 
single spin state fermions resonantly interacting in p-wave 
via a two-channel Feshbach resonance, in the sector of 
total angular momentum one. 

The central model that wc used to describe the inter- 
action depends on three parameters, the scattering vol- 
ume Vs which diverges on resonance, an effective range 
parameter a and the spatial range of the interaction b. 
Whereas b is of the order of the van der Waals length, the 
parameter a on resonance can range from a strictly posi- 
tive minimal value of the order of 1/6 up to plus infinity, 
the minimal value being model-dependent and equal to 
l/(b-!T^^'^) for our Gaussian cut-off function. In present 
experiments one estimates afe « 3 for '"'K [25] and a6 « 3 
for ^Li [30] on resonance. The two-body scattering ampli- 
tude for Vs = oo is /fe ~ — 1/a for low relative momenta 
fc ^ 1/6, so that |/fc| right on resonance is at most of the 
order of b, which is extremely small as compared to the 
usual s-wave unitary limit. As a consequence, the res- 
onant three-body problem has very different properties 
from the s-wave one. 

First, it does not exhibit the Efimov effect but it ad- 
mits two trimers, one with even parity and the other with 
odd parity, for low enough values of ab. Since the con- 
sidered sector is of angular momentum one, each trimer 
is three- fold degenerate. For Vs/b^ large and negative, 
our model Hamiltonian does not have a two-body bound 



state, so that these trimer states are examples of Bor- 
romean states. However, we estimate that the sponta- 
neous decay rate Fioss of the trimers, due to the forma- 
tion of deeply bound dimers present in current experi- 
ments with real atoms, eventually becomes larger than 
the binding energy of the trimers (over h) if one gets 
very close to their formation threshold: in the limit of a 
vanishing trimer binding energy, Fioss tends to a non-zero 
limit for the even trimer, and vanishes as the square root 
of the trimer binding energy for the odd trimer. 

Second, the atom-dimer scattering length a^d, char- 
acterizing the low-energy (that is s-wave) scattering of 
an atom on a dimer, assumes small values, of the order 
of b (or even below in absolute value for ultra-narrow 
Feshbach resonances) , except close to the odd dimer for- 
mation threshold where it diverges. The fact that Oad 
depends on the interaction range b shows that it is not 
a 'universal' quantity and it is sensitive to the micro- 
scopic details of the interaction. Furthermore, the loss 
constant in the inelastic atom-dimer s-wave scattering 
(due to the formation of deeply bound dimers) is propor- 
tional to hb/m, away from the trimer threshold, so that 
the inelastic rate may dominate over the elastic one. A 
similar conclusion was reached for the elastic vs inelastic 
scattering of p-wave weakly bound dimers [44] . We have 
also studied the atom-dimer scattering when the incom- 
ing relative wave is a p-wave: in the considered sector of 
total angular momentum one, the corresponding atom- 
dimer scattering volume V^'^ is shown analytically to be- 
come proportional to the atom-atom scattering volume 
Vs away from the even trimer formation threshold, see 
(66). This asymptotic result looks 'universal', since it 
does not involve the interaction range b. 

The recombination rate of three atoms into weakly 
bound dimers, calculated in this work in the limit of low 
relative atomic wavevectors k <^ {aVs)~^^^, has prop- 
erties more similar to the s-wave case. What remains 
specific to the p-wave case is that the rate is propor- 
tional to the square of the mean kinetic energy per par- 
ticle, see (81), as already pointed out in [42]. Apart 
from that, it includes as a factor the recombination con- 
stant /Crec, which is large close to the resonance: it di- 
verges as Vs^^/a^^^ in the large scattering volume limit, 
see (90), an asymptotic expression valid away from the 
even trimer formation threshold and which is 'universal', 
since it does not involve the interaction range b and it 
is not sensitive to the choice of the cut-off function x(k) 
of the two-channel model. In the large scattering vol- 
ume limit, the recombination constant towards weakly 
bound dimers dominates over the recombination constant 
towards deeply bound dimers, which scales as only 
(still in the limit of low relative atomic momenta with 
respect to l/(a|Vs|)^/^). If one applies this last result to 
a degenerate macroscopic gas, with a Fermi momentum 
kp w 1/(q;|Vs|)^^^, one finds a number of recombination 
events to deeply bound dimers per unit of time and vol- 
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umc scaling as 



O ( -b^n^ 



(150) 



where we used 1/a = 0(b) and n is the gas density. 

In the last part of the paper, we have made the mod- 
elization more realistic by including a fourth parameter, 
a direct attractive interaction between atoms in the open 
channel. Physically, this interaction is supposed to be 
not resonant so that it has a weak background scattering 
volume V^8, of the order of b^ and much smaller than V^. 
To stay in the regime where no deeply bound dimers ex- 
ist in the Hamiltonian on resonance, one further imposes 
V^^ < 0. We then find that the existence of the trimers 
is preserved. They remain weakly bound in the vicinity 
of some threshold values of a, provided that IV^^I does 
not exceed a few b^ ; these threshold values for a however 
depend on b and are not 'universal'. The atom-dimer 
scattering length a^d is significantly changed by the open 
channel interaction in the limit of ultra-narrow Feshbach 
resonances, where it now tends to zero for large ab as 1/a, 
rather than as l/{a^b) in our 3-parameter model; this 
confirms the non-universal character of Oad- On the con- 
trary, in presence of open channel interactions, the same 
asymptotic expressions (66) for the atom-dimer scatter- 
ing volume VJ"* and (90) for the recombination constant 
/Crec to weakly bound dimers are obtained as in our 3- 
parameter model, in terms of Vs, and in terms of Vs and 
a respectively, which confirms the 'universal' character 
of these results. 
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APPENDIX A: INTEGRAL EQUATIONS FOR 
THE A'-DEPENDENT PART OF THE ANSATZ 
FOR (3 



When injecting the ansatz Eq.(41) or Eq.(42) in the ho- 
mogeneous part of the equation Eq.(40) for /3, one faces 
the calculation of the following angular averages over the 



direction of k, 
Io(K,fc) = 

X 

Ii(K,fc) = 



^ I K + k 



K-e^)-fi(k-e^) ^_K.kb^ 



— ^ -K + k 

47r \ 2 



(Al) 



K 



-K 



l2(K,fc) 



g2 + 


A'2 + /jS + K • k 








(iK..) 


7 4^ 





(A2) 



X (k • e^) 



1 I K k 

2 'T" fc^ 



,-Kk6; 



g2 + if 2 _^ fc2 + K • k 



A3) 



Note that, for the calculations performed in this paper, 
one can restrict to the case of a non-positive total energy 
E so that we have set E = —Ti^q^/ra, q> 0, and we have 
omitted iO+. 

To perform this angular integration, we use spherical 
coordinates of polar axis K/AT. We need to evaluate first 
the integral over the azimuthal angle if, 

Bi = / ^k = fccosS^ (A4) 
Jo 27r A 

B2= J(k.e„)k = Vfc2s^^.(E, -ejE, (A5) 

where a stands for y or z, {£^=1^2.3} is an orthonormal 
basis with E3 = K/AT and 



Jo 27r 



(A6) 



[(1 - cos^ 0)% + (3cos2 9 - 1)5^36,3] (A7) 



This leads to 



Bo 



fc2 

2 



(l-cos2 6i)e„ + (3cos2 6i-l 



(K ■ e^)K ' 



(A8) 

The integration over the polar angle 9 involves basi- 
cally the integral 



^ dx x"e *^ 



2 V + X 

with V > 1, t > and integer n. The result is 



Co = —[Eiivt-t)-E,{vt + t)] 



Ci = -vCo+jo{it) 
C-2 = -vCi^iji{it) 



C3 



-VC2 - j2{it) - -ji{it) 



(A9) 



(AlO) 

(All) 
(A12) 

(A13) 
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where El (z) is the exponential integral 

n+oo „ — sz 



with 



ds ■ 



(A14) 



and jn(z) are the usual spherical Bessel functions. A 
straightforward integration over the 9 angle leads to 



Io(K,fc) 
Ii(K,fc) 

l2(K,fc) 



+ 



-(Co-C2)e, 

|(3C,-Co) + ^Ci + i-Co 



(K • e,)K 

]f2 



:(C2-Co) 



K 



K e„ 



K 



K 



+ 



n 1 

^(C^0-C2) + -(Ci-C3) 



J(3C2 - Co) + i(3C3 - Ci) + \ci + i^C2 
(K • e,)K 



k 
K 



K' 



Kk 



t = h^Kk (A18) 



X2 



{A15) 
(A16) 

For the odd sector, £' < so the source term Aq van- 
ishes. Projecting Eq.(40) onto the two components 
and (K • ezYK/ K'^, we obtain an integral system of cou- 
pled equations for Bl=o{K) and Bl^2{K), given in the 
main text, see Eq.(46), where we have introduced the two 
(A17) by two matrix 



J 



M{K, k) 



k 

' AK 



4^ (Co - C2) 



(3C2 - Co) 



fCi 



K k (or^ 

'^0 TK\'i^2 



"ifeICo - C2) - \{Ci - C3; 



2k 



Co) 



ICa 



K r< 

21^2 



(A19) 



r 



The resulting integral operator can be made hermitian 
by the change of variables 



with the two by two transformation matrix 

P : 



21/4 
_2-l/4 2-1/4 



This results in the hermitian integral equation 

4 r+oo 



(A20) 



(A21) 



= D{K) 



boiK) 
h2{K) 



dk 



60 (fc) 
&2(fc) 



(A22) 



with the two by two matrix N{K,k) = PM{K,k)p-^ 
satisfying N^{K,k) = N{k,K). For the even sector, the 
integral equation for i?L=i is given directly in the main 
text, see Eq.(76). 

To conclude this appendix, we briefly explain how to 
normalize the state vector of the trimer (when it exists) . 
The normalization can be done directly in momentum 
space by integration over internal variables, that is af- 
ter having singled out the total momentum variables Q. 



In the sector of (35) with one molecule, using the fact 
that the parameterization of the molecular and atomic 
momenta k,noi = Q/2 + K, kat = Q/2 — K, has a unit 
Jacobian, we find 



d^K 

(2^ 



l^(K)|^ 



(A23) 



In the purely atomic sector of (35), using the fact that 
the parameterization of the atomic momenta ki = Q/3-I- 
K/2 + k,k2 = Q/3 + K/2 - k, kg = Q/3 - K, has a 
unit Jacobian, using Wick's theorem and the fact that 
A(K,k) is an odd function of k, we obtain 



A*(K,k) 



2A* 



(A24) 



Using the specific form of the ansatz (42) and (41) for 
(3, with Bl=o{K), Bl=i{K) and Bl=2[K) real, and the 
link (38) between A and (3 written here for Aq = 0, one 
can get integrals of lower dimensions: 
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\ inol I mol / 



dR 



2B{K)B{k){Do - D2){K, fc)e-i(^'+'=')^' 



Gn'^h^ Jo 

dK 
3\/27r2 



6V2nHi^ Jo 



|fc4b2(A')e- 



MK) ■ Nd{K, fc)b(fc)e-i(^'+'=')^' 



(A25) 
(A26) 
(A27) 
. (A28) 



We have set E 



Da = 



n\l,Jm, B{K) = 
^ dec x"e~*'^ 



KBl=i{K) and 
(A29) 



which is minus the derivative of C„ with respect to v for 
fixed t and obeys the recursive relation Dn+i ~ C„— wZ?„. 
Also, the vector h{K) has components ho{K),b2[K), and 
the two by two matrix Nd{K, k) is obtained in replacing 
each C„ by £)„ in the matrix N{K, k). 



APPENDIX B: PRESCRIPTION FOR THE 
RECOMBINATION RATE OF FERMIONS 

We wish to derive the formula (81) giving the rate of 
dimer formation in a gas of fermions at low kinetic en- 
ergy, in terms of a recombination constant depending on 
the interaction and the expectation value (. . .)o of some 
operator in the unperturbed state of the gas. We start 
with the intuitive idea that a dimer formation can take 
place by three-body collision when the mutual distances 
of the atoms are at most of the order of the dimer radius 
a, hence the heuristic formula 



— iVdim (X 

at 



d^ri d^r2 d^r^ e 



x(^t(ri)Vnr2)V'Hr3)V'(r3)^/'(r2)V'(ri))o. (Bl) 

We have taken for simplicity a Gaussian cut-off function, 
in terms of the hyperradius p defined as 



= ^(r,-R)2 = i(ri-r2)^ 
1=1 



^[r3-(ri+r2)/2]^ 
(B2) 

where R — X]i=i r^/S is the center of mass position. Un- 
der the assumption that the typical wavcvectors popu- 
lated in the uncorrelated state of the gas are much smaller 
than 1/cr, which is the case here since we assume rela- 
tive momenta < (jdim, we can expand each field operator 
^^{vi) around the center of mass position R of the three 
Yi in powers of — R. Since ?/'2 = one has to go to 
second order in 5r: 



■0(r3)'0(r2)-0(ri) 



^(R) 



where a and /3 run over the three directions of space x, 
y and z, and the matrix A is given by 



(B4) 



We restrict to this leading order in 6r. It remains to 
integrate the Gaussian weighted products AapA^s over 
the internal variables for fixed center of mass position R, 
by using the Jacobi coordinates. Since Aaf3 for a ^ (3 
is odd with respect to the reflection along direction a or 
along direction /3, this integral vanishes if {a, (3} ^ {7, 5}. 
The invariance of the integral by permutation of the x, 
y and z axis leads to the final prescription (81). 



APPENDIX C: RECOMBINATION CONSTANT 
CLOSE TO THE TRIMER FORMATION 
THRESHOLD 



The goal is to derive the approximate formula (93) giv- 
ing the recombination constant /Crcc for large scattering 
volumes and for a value of arcs close to the threshold for 
the even trimer formation. 

To this end, we rewrite the integral equation (76) 
as a sum of its = 00 value and a remainder, then 
we treat the remainder perturbatively. Taking as un- 
known F{K) = Kg{K)/Vs, using Eq.(26) and the iden- 
tity l/{X-iO+) = P^+iTrS{X), where 7^ is the Cauchy 
principal value, we obtain the rewriting 

-^F{K) ™ Io[F] (K) - h [F] (K) ~ tA{K)F{K) = S{K) . 



We have introduced the two functions 



(CI) 



AiK) = if(Co-C2)(ir,i^dim)e-^'''(^'+^-»)(C2) 



SiK) 



mA 

and the two operators 



Kok^e-'"''''/' 



(C3) 



(B3) 
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h{K] 



F{K) + - I dK' §j{Co - C2){K,Ky~l'''^^'+^''^F(K% 



K 



(C4) 



A'2 - K'i 



F{K) + r°°dK' ^{Co - C2){K,K')e-i'''(^'-+'<''^V 

TT Jo K' 



F{K' 



X'2 



^C5) 



r 



(C6) 



The background part of the solution. F\^g{K), and its 
derivative F^^(K), are supposed to be uniformly bomided 
in Vs and ares, even in the vicinity of ofj-cs = Qiti™"- From 
the low K behavior of Co{K, K') - C2{K, K'), one finds 
that the solution F{K) satisfies F(0) = 45(0)/aros, this 
value being reached quadratically in K. Since \F'\ is 
not uniformly bounded in Vs and Ores, this does not give 
information on the value F(i^dim)- One has also -Fbg(O) = 
4S'(0)/aros, but for a |F^g| bounded by Cbg, where the 
constant Cbg does not depend on Vs or arcs, one now has 



|J^bg(if, 



dim ) 



< 



Cbg-ft^dir 



Fhg{Kdim) 



where we have set h{K) ~ qcxp{q^b^) crk{qb) with q = 
V3K/2. 

The operator J such that J[G]{K) = h[F][K)/K 
where F{K) = KG{K), is self-adjoint and thus has real 
eigenvalues; wc find numerically that it has one and only 
one discrete eigenvalue, that we called q;5;™"/4, with the 
corresponding normalized eigenvector K — > go{K). As a 
consequence, the operator Iq admits, with the eigenvalue 
a°™"/4, a discrete eigenvector uo{K) = Kgo{K), with F{Kdn-a) 
the corresponding adjoint (left-eigenvector) vo{K) = 
go{K)/K. We note that the kernel in /q behaves as 
K'"^ /{K"^ + K'^) at low momenta, which is bounded, so 
that vq{K) is bounded, go{K) vanishes linearly with K in 
K = and uq{K) vanishes quadratically [78]. Apart for 
this discrete eigenvalue, we numerically find that Iq has 
a continuous spectrum extending from to 1/ {Ay/nb) . 

For Qfios close to a"™", in the large scattering vol- 
ume limit, a small denominator (of the order oi 6 = 
ares — c^th"^) appears in the direction of uq when one 
solves (CI). This small denominator is weakly perturbed 
by Ii (which shifts the value of arcs for which the de- 
nominator has a minimal modulus) and by the imaginary 
part involving the function A(K) (which prevents the de- 
nominator from exactly vanishing). These effects can be 
included systematically by using the ansatz 



where we used (i'o|/o|-Fbg) = (a^J"V4)(t;o|Fbg) (after 
justification). Here Dirac's notation means {f\g) = 

J^°° dK f* {K)g{K). Replacing cq by its expression in 
(C6); setting K — Kdim and putting all terms on a com- 
mon denominator, we see that the imaginary contribu- 
tion iuo(Kdim){'Vo\A) exactly cancels in the numerator. 
We finally obtain the still exact expression 



= T[S/4-{vo\Ii\uo)-tuo{Kdin.){vo\A)r^ 
X {S/4-{vo\h\uo)/T + uo{Kdi^-,)[{vo\S) 

+ («o|/l|^bg)]/[Fbg(A'dim)-^]}, (C9) 



where = 1 - (uo|-Fbg)uo(-Adim)/-Fbg(/idim)- 

The last step is to expand the various terms to leading 
order in if dim- Expanding Ii to leading order in i^dim and 
using the fact that uo is an eigenvector of Iq to simplify 
the integral expression appearing in this leading order 
form of Ii , we obtain 



lim 



{vo\I: 



i\uo) 



"dim 



{vo\vo) > 0, 



(CIO) 



and the quantity {vq\vo) is readily evaluated numerically. 
This gives a position of the peak in /Croc shifted to a value 
of arcs larger than a^^°'^ by a 0{K^^^^), see (95). Since 
uq{K) vanishes quadratically in A', one has wo(A'dim) 
of the order of K\^^^; since A{K) vanishes linearly with 
ATdim for a fixed AT, one has (wo|^) of the order of ATdim- 
A numerical calculation of the corresponding coefficients 
leads to (96). Amusingly, using the low- A' expansion of 
Co (AT, AT') — C2{K, AT'), we find the mathematical cquiv- 



(Cll) 



(C7) which leads to (97). Another result is 



This is the property of Abg that we shall need. 

Injecting the ansatz (C6) into (CI) and projecting onto 
Uq by using the adjoint vector vq, we obtain the exact 
expression 



[S/A - (uo|/i|uo) " iuo{Kdi, 
-(V4)(«o|Abg) + (i'o|/i|Abg) 



,)MA)]co = {vo\S) 
-iAbg(A-dim)(«o|^),(C8) 



(^^ol/ilAbg) -0(ifd'imlnA'd 



(C12) 



so that this contribution, being multiplied by wo(Ardim) 
in (C9), may be neglected at this order. Finally, we note 
that Abg (ATdim) and differ from 45(0)/arcs and 1 re- 
spectively by terms of order A'd 
neglect to obtain (93). 



and AT^inj, that we 
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